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ABSTRACT 


The  eclipsing  spectroscopic  binary  system  DE 
Draconis  is  re-examined,  using  86  new  observations  ob¬ 
tained  between  1976  and  1978  on  the  1.8m  telescope  at  the 
Dominion  Astrophysical  Observatory  in  Victoria  and  one 
new  point  from  the  David  Dunlap  Observatory  in  Toronto, 
together  with  18  older  data  from  the  DDO  in  Toronto.  The 
52  observations  outside  eclipse  allow  a  more  reliable  de¬ 
termination  of  the  orbital  elements  to  be  made,  using  the 
Lehmann-Filhes  method  of  solution,  and  this  allows  a 
velocity  curve  for  the  primary  to  be  drawn.  The  spectro- 
graphic  results  are  combined  with  some  previously  pub¬ 
lished  photometric  and  four-colour  results  to  give  an 
approximate  picture  of  the  physical  parameters  of  this 
system. 

The  very  detailed  spectrographic  study  of  the 
velocity  changes  of  the  primary  star  on  53  plates  taken 
during  eclipse  clearly  shows  the  existence  of  a  rotation 
effect  in  the  primary's  velocity  curve.  This  is  analysed 
using  two  methods,  by  Petrie  and  Hosokawa.  The  rotational 
velocity  of  the  primary  star,  spectral  type  B9V ,  is 
found  to  be  about  180  km/s,  indicating  that  the  star  does 
not  rotate  in  synchronism  with  the  system.  All  results 
support  the  previous  conclusion  that  the  secondary, 
although  not  directly  observed,  is  a  main  sequence  G-type 
star  with  normal  characteristics. 
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CHAPTER  I 


INTRODUCTION 

Historical  Recognition  of  Binary  Stars 

The  term  'binary  star*  was  first  used  by  W.  Her- 
schel  in  1802  to  describe  "...  a  real  double  star — the 
union  of  two  stars,  that  are  formed  together  in  one  sys¬ 
tem,  by  the  lav/s  of  [gravitational]  attraction  ..." 
(Herschel  1802).  However,  the  term  'double  star'  had  been 
used  much  earlier,  by  Ptolemy  in  describing  v  Sagitarii — 
two  fifth  magnitude  stars  with  an  angular  separation  of 
about  14  minutes  of  arc  (somewhat  less  than  half  the  ap¬ 
parent  diameter  of  the  moon) . 

Thus  not  all  'double  stars'  are  binary  systems. 

The  first  binary,  discovered  in  1650,  was  £  Ursae  Majoris 
(Mizar) .  It  was  also  the  first  double  star  to  be  observed 
photographically,  in  1857,  and  its  principal  component 
was  the  first  spectroscopic  binary  to  be  identified,  in 
1889.  Algol  (3  Persei)  became  the  second  known  spectro¬ 
scopic  binary  and  the  first  of  the  eclipsing  binaries, 
when  its  spectrum  was  investigated  because  of  the  periodic 
loss  of  light  observed  by  Goodricke  (1783) . 

The  earliest  discoveries  made  in  the  seventeenth 
and  most  of  the  eighteenth  centuries  were  all  accidental, 
with  no  thought  that  double  stars  were  more  than  optical 
coincidences.  Halley's  observation  in  1718  that  some 


1 


. 


2 


stars  were  actually  in  motion  probably  led  to  speculations 
in  the  middle  of  the  eighteenth  century  on  the  existence 
of  systems  of  stars.  In  1767  John  Mitchell  pointed  out 
that  the  frequency  of  occurence  of  double  stars  was  much 
too  high  to  be  accounted  for  by  the  probability  of  chance. 
Mitchell's  younger  contemporary,  W.  Herschel,  remained 
unconvinced,  but  finally  in  1803  (Herschel  1803)  conceded 
that  some  double  stars  were  indeed  true  binaries. 

After  Herschel,  the  list  of  visual  double  stars 
continued  to  grow  through  observations  made  by  many 
astronomers,  notably  F.G.W.  Struve  in  Russia  and  S.W. 
Burnham  in  the  United  States.  These  wide  pairs  of  stars 
(i.e.  resolvable  through  a  telescope)  for  which  orbital 
elements  and  trigonometric  parallaxes  are  available,  are 
important  sources  for  the  determination  of  stellar  masses. 
These  stars  do  not  influence  each  other  except  for  their 
mutual  gravitational  attraction,  and  thus  behave  as  mass 
points.  Of  increasing  interest  to  the  modern  astronomer 
are  those  pairs  in  which  the  separation  of  the  components 
becomes  roughly  comparable  to  the  stellar  size — these  are 
termed  'close'  pairs,  and  cannot  be  visually  resolved. 
However,  other  methods  of  observation  have  revealed  many 
of  these  close  binaries. 

In  the  latter  part  of  the  seventeenth  century  the 
star  Algol  (A1  Ghul,  as  it  was  named  by  the  people  of  the 
ancient  Near  East)  was  noted  as  having  variable  brightness. 


3 


Finally  in  1782-3,  John  Goodricke  observed  Algol,  realized 
the  short  periodicity  of  its  brightness,  and  postulated 
that  this  variability  was  due  to  "...  the  interposition 
of  a  large  body  revolving  around  Algol..."  (Goodricke  1783). 
However,  only  in  1889  (Vogel  1890)  when  Algol  was  recog¬ 
nized  as  a  spectroscopic  binary  having  times  of  conjunction 
coincident  with  times  of  minimum  light,  were  Algol  and 
similar  variable  stars  recognized  as  binary  systems. 

3  Lyrae  was  also  discovered  by  Goodricke  in  1784  to  be  a 
light  variable,  and  was  thus  the  second  eclipsing  binary 
to  be  recognized. 

Classification  of  Binary  Stars 

There  are  three  classes  of  binary  stars,  and  a 
given  system  may  be  a  member  of  one,  two,  or  even  three 
of  the  classes.  The  first  class  is  composed  of  systems 
whose  components  are  far  enough  apart  to  be  optically 
resolved,  called  ’visual'  binaries.  Direct  observations 
of  their  apparent  positions  can  be  made  with  a  telescope 
and/or  on  photographic  plates.  Other  double  stars  are 
either  too  far  away  from  Earth  and/or  too  close  to  each 
other  to  be  visually  double.  The  second  class  can  be  de¬ 
tected  by  the  wavelength  shifts  in  spectra  giving  a  per¬ 
iodic  variation  of  the  radial  velocity  (or  velocities) , 
and  are  called  'spectroscopic'  systems.  The  third  class 
is  identified  by  photometric  observations  of  periodic 
variations  in  the  light  emitted,  and  are  called 
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'eclipsing'  binaries.  These  three  classifications, 
however,  are  not  mutually  exclusive.  For  example,  a 
spectroscopic  binary  which  has  an  appropriate  orbital  in¬ 
clination  can  in  principle  also  be  observed  as  an  eclip¬ 
sing  system.  The  type  assigned  a  binary  system  in  this 
way  is  thus  dependent  on  its  distance  from  Earth,  and 
there  is  much  overlap.  However,  any  binary  system  is 
completely  specified  by  seven  parameters:  the  mass, 
radius,  and  luminosity  of  each  component,  and  the  mean 
distance  between  their  minimum  and  maximum  separation. 

A  different  method  of  classification  was  proposed 
by  Kopal  (1955)  ,  based  on  the  theoretical  upper  limit  to 
the  size  a  star  can  reach  in  the  presence  of  a  companion. 
This  size  is  usually  thought  of  as  the  first  of  the  zero- 
velocity  curves  of  the  restricted  three-body  problem, 
called  the  critical  surface,  and  the  maximum  volumes  at¬ 
tainable  are  the  'Roche  lobes'.  The  relative  dimensions 
are  completely  specified  by  the  ratio  of  the  two  masses, 
and  systems  are  classified  according  to  whether  or  not 
one  or  both  of  the  components  fill  their  respective  Roche 
lobes.  This  classification  uses  information  on  the 
masses  and  radii  of  the  components. 

The  term  'close  binary  system'  (as  opposed  to 
'wide')  usually  refers  to  one  in  which  one  component 
affects  the  evolution  of  the  other,  and  thus  is  closely 
related  to  Kopal' s  classification.  In  normal  stellar 
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evolution,  the  radius  of  the  evolving  star  increases,  and 
thus  in  a  binary  system  the  natural  radius  may  attempt  to 
exceed  its  Roche  lobe,  so  that  the  two  stars  will  affect 
each  other's  evolution.  In  this  sense  close  binaries  en¬ 
compass  visual,  spectroscopic,  and  eclipsing  systems,  and 
are  directly  involved  with  the  evolutionary  history  of 
binaries.  The  variety  of  observed  binaries  is  due  to  the 
fact  that  at  different  stages  one  of  the  components  (the 
more  massive)  fills  its  Roche  lobe. 

The  recognition  of  a  visual  binary  depends  on  the 
separation  of  its  components  and  its  distance  from  the 
Sun.  Thus,  visual  binaries  must  be  relatively  near  neigh¬ 
bours,  within  a  few  hundred  parsecs.  Spectroscopic 
binaries  can  be  recognized  out  to  a  few  thousand  parsecs 
with  high  dispersion  spectrographs.  Beyond  this  limit 
only  eclipsing  binaries  can  be  detected  (if  the  eclipses 
are  fairly  deep)  and  thus  their  recognition  in  galactic 
clusters  can  aid  in  the  determination  of  the  distance  to 
the  particular  cluster.  If  the  age  of  the  cluster  can  be 
found  from  its  colour-magnitude  diagram,  a  check  can  be 
made  on  theories  of  evolution  of  close  binaries.  Spec¬ 
troscopic  binaries  which,  due  to  orbital  inclination,  are 
also  eclipsing  have  been  the  main  source  of  information 
on  individual  stellar  dimensions  and  masses,  leading  to 
the  empirical  mass-luminosity  relationships.  If  the 
parallax  of  a  system  is  also  known,  the  effective  temper- 
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atures  of  the  stars  can  be  determined.  Such  data  provide 
the  basis  for  studies  of  stellar  structure. 

Orbital  Elements 

The  orbit  of  a  binary  system  is  completely  defined 
by  seven  elements,  of  which  the  period  is  the  most  impor¬ 
tant  because  the  generalized  form  of  Kepler's  Third  Law 
relates  the  period  P  to  the  semi-major  axis  a,  and  to  the 
stellar  masses: 

3 

~  «  (mi+  M2)  (1.1) 

Thus,  the  following  quantities  will  completely  specify 
the  system  (Batten  1973): 

P  :  the  orbital  period,  expressed  in  days  for  spec¬ 
troscopic  or  eclipsing  binaries,  or  in  years  for 
visual  binaries 

i  :  the  inclination  of  the  orbital  plane  to  the  tan¬ 
gent  plane  to  the  sky 

Q,  :  (applicable  only  to  visual  binaries)  the  position 
angle,  measured  from  north  towards  east,  of  the 
line  of  nodes  joining  the  intersection  of  the  or¬ 
bital  and  tangent  planes,  measured  in  the  tangent 
plane 

a)  :  the  longitude  of  periastron,  which  is  the  angle 

between  the  directions  of  the  ascending  node  (where 
the  star  crosses  the  tangent  plane  moving  away  from 
the  observer) ,  and  of  the  point  of  periastron 
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(the  point  of  closest  approach  of  the  two  stars) . 

It  is  measured  in  the  orbital  plane  in  the  direc¬ 
tion  of  orbital  motion.  Visual  binaries  have  oo 
given  for  the  secondary  (fainter)  component,  but 
spectroscopic  and  eclipsing  systems  specify  oo  for 
the  orbit  of  the  primary  component — the  two  differ 
by  180° 

a  :  the  semi-major  axis  of  the  relative  orbit  (usually 
in  astronomical  units  or  kilometers) 
e  :  the  eccentricity  of  the  orbit,  ranging  between 

zero  (for  a  circular  orbit)  and  one  (for  a  straight 
line) 

T  :  the  time  of  periastron  passage,  usually  expressed 
as  a  Julian  date 


Fig.  1.  Orbital  elements  (from  Batten  1973,  p.  8) 
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Determination  and  Significance  of  Orbital  Elements 

For  a  visual  binary  the  elements  P,  i,  u,  e,  and 
T  can  be  found.  The  ascending  and  descending  nodes  can¬ 
not  be  distinguished  unless  there  are  also  radial  velocity 
measurements,  so  there  is  an  uncertainty  of  180°  in  ft. 
Usually  i  is  taken  to  be  between  0°  and  90°  if  the  appar¬ 
ent  motion  is  direct  (i.e.  in  the  direction  of  increasing 
position  angles);  otherwise  i  lies  between  90°  and  180°. 

The  ascending  node  is  assumed  to  be  that  node  for  which 
ft  <  180°.  The  apparent  semi-major  axis  in  seconds  of  arc 
can  be  derived,  but  an  absolute  value  is  not  available 
unless  the  distance  to  the  system  can  be  independently 
determined.  Radial  velocity  data  can  give  the  scale 
(relative  size)  of  the  orbit,  and  thus  the  distance  can 
be  determined.  Using  Kepler's  Third  Law,  assuming  normal 
masses  for  the  spectral  type  of  each  component,  the  major 
axis  (and  the  distance)  can  be  estimated.  This  distance 
estimate  expressed  as  a  parallax  is  called  the  dynamical 
parallax.  The  periods  of  visual  binaries  are  usually 
very  long  (years) ,  and  so  period  determinations  are  often 
uncertain,  based  on  only  one  arc  of  the  curve. 

In  addition,  for  a  visual  binary,  the  relative 
brightness  of,  or  the  difference  in  magnitude  between, 
widely  separated  components  can  be  measured  photometri¬ 
cally  but,  until  recently,  the  systems  with  orbital  periods 
short  enough  to  allow  other  element  determinations  had 
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such  small  angular  separation  that  accurate  magnitude  dif¬ 
ferences  could  not  be  obtained. 

Elements  determined  for  spectroscopic  binaries 
are  different,  and  the  amount  of  information  depends  on 
whether  the  spectra  of  both  components  are  visible  (a 
double-lined  spectroscopic  binary) ,  or  the  spectrum  of 
only  the  brighter  component  (a  single-lined  spectroscopic 
binary)  is  visible.  Usually  if  the  two  stars  differ  by 
more  than  about  one  magnitude,  only  the  spectrum  of  the 
brighter  star  is  visible.  The  period  can  be  determined 
(see  Chapter  3)  for  either  case,  and  from  the  shape  of 
the  velocity  curve  the  eccentricity,  e,  and  the  longitude 
of  periastron,  oo,  can  be  found.  In  addition,  y,  the 
radial  velocity  of  the  centre  of  mass  of  the  system  (also 
known  as  the  systemic  velocity) ,  and  K^,  half  the  total 
range  of  the  radial  velocity  variation  of  the  brighter 
companion,  are  obtained  from  the  velocity  curve.  The 
radial  velocity  reaches  its  maximum  at  the  ascending  node 
and  its  minimum  at  the  descending  node.  The  observed 
velocity  contains  the  velocity  of  the  system  as  a  whole 
as  well  as  the  velocity  of  the  individual  components. 

From  the  definition  of  ,  the  quantity  a^sin(i)  can  be 
calculated,  where  a^  is  the  semi-major  axis  of  the  abso¬ 
lute  orbit  of  the  brighter  star.  In  addition,  if  the 
spectrum  is  double-lined,  K2  and  a2sin(i)  for  the  fainter 
star  can  also  be  found.  Using  Kepler's  Third  Law,  for  a 


' 

. 


single  spectrum  binary,  a  quantity  called  the  'mass  func¬ 
tion',  f (M) ,  can  be  calculated: 
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f  (M) 


(m1+m2)  2 


.  3 
sm 


(i) 


1.035x10  Vp(l-e2)3/2  M  (1.2) 
1  © 


where  P  is  in  days  and  K  is  in  km/s. 
will  give  the  ratios 


a2  k2  M]_ 


Double-lined  spectra 


(1.3) 


and  thus  the  reduced  masses 

M,  0sin3(i)  =  1.035x10_7Ko  .  (K. +K0 ) 2p (1-e2 ) 3/2  M  (1.4) 

can  be  determined.  Here  (K  +K2 )  is  proportional  to  the 
semi-major  axis  projected  onto  a  plane  containing  the 
line-of-sight ,  and  (a^+a  )=a.  If,  as  is  generally  the 
case,  the  fainter  star  is  also  the  less  massive,  then 
minimum  values  of  the  reduced  masses  are: 


3 

M1  2sin  (i)  =  4f]_  2  ^ 


(1.5) 


For  double-lined  spectroscopic  binaries,  the  rel¬ 
ative  intensities  of  the  lines  are  directly  proportional 
to  the  stellar  luminosities.  Petrie  (1939;  1950)  developed 
methods  of  measuring  the  intensities  and  thus  obtaining 
the  luminosities  via  the  relative  depths  or  equivalent 
widths  of  selected  spectral  lines.  However,  difficulties 
may  arise  if  the  two  stars  are  of  very  different  spectral 
types,  or  if  the  line  profiles  in  the  spectra  are  of  dif- 
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ferent  shapes,  or  if  (as  in  early-type  stars)  the  lines 
are  very  broad  and  so  not  properly  resolved. 

From  the  light  curves,  the  periods  can  be  found 
for  eclipsing  binaries.  The  minimum  distance  between  the 
two  stars  projected  onto  the  plane  of  the  sky  is  related 
to  the  inclination  of  the  orbit,  which  affects  the  shape 
and  depth  of  the  light  curve,  The  angle  0,  is  indetermin¬ 
ate,  since  the  direction  of  motion  of  the  eclipsing  star 
cannot  be  found  (except  through  radial  velocity  measure¬ 
ments)  .  The  absolute  size  of  the  orbit,  a,  cannot  be 
found  unless  spectroscopic  data  are  also  available.  How¬ 
ever,  the  shape  and  length  of  the  eclipses  depend  on  the 
relative  radii  of  the  stars  in  terms  of  their  separation, 
and  these  fractional  radii  can  be  determined  from  the 
light  curve.  If  spectroscopic  observations  are  available, 
the  absolute  radii  can  be  calculated.  If  both  eclipses 
are  observable,  the  displacement  between  them  as  a  frac¬ 
tion  of  the  period  is  proportional  to  ecos  (w)  to  a  first 
approximation  when  e  is  small.  Short  period  eclipsing 
binaries  usually  have  nearly  circular  orbits.  For  circu¬ 
lar  orbits  or  elliptical  orbits  with  the  major  axis 
along  the  line-of-sight ,  ecos  (oo)  =  0,  and  the  secondary 
eclipse  occurs  half-way  between  two  primary  eclipses. 

The  difference  in  the  durations  of  the  primary  and  secon¬ 
dary  eclipses,  relative  to  the  period,  is  proportional  to 
esin(m)  when  e  is  small.  Thus,  e  and  to  can  in  principle 
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be  separately  determined,  but  usually  the  value  of  esin(rn) 
has  a  large  uncertainty. 

For  an  eclipsing  binary,  the  combined  light  of  the 
two  stars  outside  eclipse  is  taken  as  the  unit  of  light. 

If  one  eclipse  is  total  (a  deep,  flat-bottomed  minimum) , 
the  relative  luminosities  are  determined,  since  the  depth 
of  totality  gives  directly  the  magnitude  difference  of 
the  two  stars.  For  partial  eclipses  (assuming  a  nearly 
circular  orbit--as  is  usually  the  case) ,  the  same  area  of 
each  star  is  eclipsed  and  the  ratio  of  the  depths  of  the 
two  eclipses  gives  the  ratio  of  the  surface  brightnesses. 
When  the  radii  are  calculated  from  the  shape  of  the  light 
curve,  the  ratio  of  luminosities  can  then  be  determined, 
expressed  at  one  wavelength.  If  measurements  are  made  at 
different  wavelengths,  the  results  can  be  used  to  find 
the  difference  in  spectral  types,  thus  adding  to  or  con¬ 
firming  the  spectroscopic  data. 

In  addition,  the  limb-darkening  of  the  eclipsed 
star  can  be  determined  by  its  effect  on  the  shape  of  the 
light  curve.  However,  accurate  calculations  are  diffi¬ 
cult,  since  it  affects  only  the  small  portion  of  the 
curve  when  the  limbs  of  the  star  are  being  covered  or 
uncovered.  Also,  in  principle,  the  distortion  of  the 
components  in  an  eclipsing  system  can  be  determined, 
since  during  times  of  eclipse  the  ellipsoidal  star  shows 
less  surface  area  than  outside  eclipse.  However,  the 
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observed  variations  in  light  between  eclipses  also  involve 
other  complications  such  as  gravity-darkening,  tidal  dis¬ 
tortion,  reflection  effects,  and  the  presence  of  circum- 
stellar  matter.  The  point  on  a  distorted  star  closest  to 
the  observer  will  not  be  the  geometric  centre  of  the  ap¬ 
parent  disc,  so  the  centre  of  maximum  light  is  displaced 
from  the  star's  centre  of  mass.  Thus,  the  light  changes 
expected  within  and  between  minima  are  very  complicated 
and  difficult  to  interpret. 

Selection  Effects 

For  visual  binaries,  the  most  important  factor 
affecting  their  discovery  is  the  distance  from  the  Sun. 
Thus,  the  observed  frequency  of  visual  binaries  appears 
to  be  higher  among  stars  of  later  spectral  type,  because 
hot,  luminous,  early-type  stars  are  scarce  in  the  solar 
neighbourhood.  Another  consideration  is  the  magnitude 
difference  between  the  components:  a  very  luminous  star 
will  more  likely  have  a  much  fainter  companion,  while  a 
fainter  star  would  have  a  companion  of  similar  luminosity. 
Therefore,  again,  visual  binaries  are  harder  to  detect 
among  the  more  luminous  stars. 

The  larger  the  changes  in  velocity,  the  easier  it 
is  to  detect  a  spectroscopic  binary.  Thus,  it  is  easier 
to  detect  spectroscopic  binaries  among  more  massive  sys¬ 
tems  which  have  larger  velocity  variations  for  a  given 
period  and  inclination,  since  orbital  velocity  is  propor- 
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tional  to  (mass)  x  .  The  more  massive  stars  are  usually 
the  hot  luminous  ones.  Systems  of  later  spectral  type 
stars  have  generally  longer  periods,  and  since  velocity 
variation  is  inversely  proportional  to  the  period,  these 
systems  would  be  more  difficult  to  detect.  However,  the 
velocity  of  early-type  stars  whose  spectra  contain  few 
diffuse  lines  is  more  uncertain  than  for  later-type  stars 
which  exhibit  many  sharp  spectral  lines.  Thus,  more  long 
period  systems  with  low  velocity  amplitudes  will  be  ob¬ 
served  among  late  spectral  type  stars.  Because  double- 
lined  spectra  can  be  observed  only  in  systems  where  the 
luminosities  of  the  two  components  are  similar  (Am  less 
than  about  one  magnitude) ,  the  masses  and  absolute  dimen¬ 
sions  are  only  available  for  those  systems  whose  mass 
ratios  are  near  unity.  Some  systems  have  orbital  elements 
(e.g.  e,  co,  K)  which  make  detection  difficult.  Therefore, 
in  any  group,  the  number  of  observed  spectroscopic 
binaries  will  be  underestimated,  an  effect  more  important 
for  early-type  stars  with  lower  accuracy  in  their  velocity 
determinations.  Double-lined  and  single-lined  binaries 
have  different  probabilities  of  discovery,  since  double- 
lined  binaries  are  detected  from  the  spectral  line  doubling 
while  single-lined  systems  are  detected  through  velocity 
changes . 

The  probability  of  detecting  eclipsing  binaries 
depends  on  the  amount  of  time  the  light  is  less  than  max- 
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imum- — that  is,  on  the  orbital  period  and  duration  of 
eclipse  as  a  fraction  of  that  period — ,  on  the  depth  of 
the  minima,  and  on  the  apparent  brightness  of  the  system. 
The  depth  and  duration  of  the  eclipse  are  determined  in 
turn  by  the  relative  radii  and  luminosities  of  the  com¬ 
ponents.  A  system  of  two  similar  stars,  which  display 
two  equal  eclipses,  is  harder  to  detect,  because  the  light 
cannot  vary  by  more  than  about  0™75,  and  also  because  the 
eclipses  will  more  likely  be  only  partial.  The  geometric 
probability,  i.e.  the  probability  that  a  system  will  be 
inclined  so  as  to  show  an  eclipse,  is  higher  for  systems 
containing  a  giant  or  super-giant  than  those  with  main- 
sequence  stars.  But  these  eclipsing  systems  containing 
widely  disparate  components  are  generally  the  longer 
period  binaries,  and  it  has  been  shown  that  long  period 
eclipsing  binaries  have  a  low  probability  of  being  photo¬ 
metrically  detected. 

Due  to  selection  effects,  double-lined  spectro¬ 
scopic  binaries  tend  to  be  found  among  the  hotter,  more 
massive  stars.  Therefore,  our  knowledge  of  masses 
greater  than  1  derives  largely  from  spectroscopic 
eclipsing  systems.  Few  stars  with  masses  greater  than 

2  M  are  found  in  visual  binaries.  The  determination  of 
o 

masses  from  eclipsing  systems  is  subject  to  the  uncer¬ 
tainty  in  sin(i).  However  i  must  be  fairly  large  to  allow 
eclipses  to  occur,  and  sin(i)  is  not  very  sensitive  to 


16 


the  value  of  i  for  i  close  to  90°.  Thus,  the  major  un¬ 
certainty  in  masses  is  due  to  the  uncertainty  in  and 


Velocity  Curve  Distortion:  The  Rotation  Effect 

Three  types  of  disturbance  in  the  radial  velocity 
curves  of  eclipsing  spectroscopic  binaries  have  been  iden¬ 
tified,  and  they  are  frequently  discovered  from  the  fact 
that  the  value  of  the  eccentricity  from  radial  velocity 
measurements  is  not  the  same  as  that  from  photometric 
observations.  The  three  effects  are  the  rotation  effect, 
the  presence  of  gas  streams  in  the  system,  and  the  reflec¬ 
tion  effect.  Of  these,  the  rotation  effect  is  of  primary 
interest  to  this  study  because  of  its  manifestation  in 
the  radial  velocity  curve  of  DE  Draconis.  This  effect 
occurs  at  the  time  of  eclipse:  as  one  side  of  the  rotating 
star  is  eclipsed,  effective  centres  of  the  spectral  lines 
are  shifted  according  to  the  direction  in  which  the  still 
exposed  limb  is  rotating.  At  mid-eclipse,  equal  areas  of 
both  limbs  are  exposed,  so  the  effect  cancels.  After  mid¬ 
eclipse  the  lines  are  shifted  in  the  other  direction,  as 
more  of  the  one  limb  is  uncovered  while  the  other  is 
eclipsed.  As  more  of  the  disc  of  the  eclipsed  star  ap¬ 
pears,  the  effect  from  each  side  cancels.  The  lines  be¬ 
come  more  symmetrical  until  outside  eclipse  the  lines  lie 
once  more  at  the  wavelength  corresponding  to  the  line-of- 
sight  velocity  of  the  star.  Spectroscopic  evidence  has 
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shown  that  stars  in  a  binary  system  generally  rotate  in 
a  direct  sense  with  an  angular  velocity  co  which  is  equal 
to  the  Keplerian  orbital  angular  velocity,  or  sometimes 
much  larger:  in  circular  orbits,  synchronism  usually 
exists  between  the  rates  of  rotation  and  revolution,  while 
in  eccentric  orbits  the  rotation  rate  is  generally  faster 
than  the  revolution  rate.  For  a  rigid  spherical  rotator 

of  radius  R  and  equatorial  angular  velocity  co,  a  quantity 

5  V 

known  as  the  'rotation  factor'  (=  — —  )  can  be  calculated 

ooR 

as  a  function  of  phase  for  a  given  orientation  of  the 
axis  of  rotation  to  the  orbital  plane,  from  elements  of 
the  system  determined  from  spectroscopic  and  photometric 
data.  The  quantity  6V  is  the  difference  at  that  phase 
between  the  observed  radial  velocity  and  the  calculated, 
or  expected,  radial  velocity  for  purely  orbital  motion. 

The  rotation  factor  is  symmetrical  with  respect  to  the 
moments  of  conjunction  when  the  eclipse  starts  and  ends, 
if  the  rotation  axis  is  perpendicular  to  the  orbital  plane. 
For  partial  eclipses,  it  vanishes  at  mid-minimum,  and  is 
equal  but  of  opposite  sign  before  and  after.  For  total 
eclipses,  it  becomes  ±sin(i)  at  the  inner  contacts.  The 
symmetry  disappears  if  the  axis  of  rotation  in  inclined 
from  the  perpendicular,  and  the  asymmetry  should  vanish 
and  reappear  periodically  in  the  opposite  direction  as 
the  axis  of  rotation  precesses  under  the  gravitational 
influence  of  the  companion. 
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From  the  numerical  calculations  of  the  rotation 
factors  and  the  observed  6V  for  each  phase,  the  equator¬ 
ial  velocity  ooR  for  the  star  can  be  calculated.  If  co  is 
assumed  to  be  the  same  as  for  the  Keplerian  orbital  motion 
(synchronous  rotation  and  revolution) ,  the  absolute  value 
of  the  radius  R  can  be  calculated.  The  axial  rotation  of 
a  star  will  cause  the  spectral  line  profiles  to  be 
broadened  by  an  amount  determined  by  the  maximum,  differ¬ 
ence  in  radial  velocity  at  opposite  limbs  of  the  star. 

If  the  stars  are  spherical,  the  broadening  is  symmetrical 
with  respect  to  the  centre-of-mass  velocity  position. 
However,  in  a  close  binary  system,  the  distortion  of  the 
stars  due  to  various  effects  causes  the  line  profiles  to 
broaden  and  become  asymmetric;  and  if  a  spherical  rotating 
star  undergoes  an  eclipse,  the  lines  will  be  asymmetric 
as  well. 


Previous  Work  on  the  Rotation  Effect 

The  fact  that  the  rotation  of  an  eclipsing  binary 
could  produce  a  measurable  effect  was  recognized  by  Holt 
(1893) ,  Schlesinger  (1909) ,  and  Forbes  (1911) .  Its  effect 
in  producing  contradictory  results  from  spectroscopically 
versus  photometrically  derived  orbital  elements  was  shown 
by  Hellerich  (1922;  1924).  He  demonstrated  that  the  dis¬ 
crepancy  could  often  be  removed  by  excluding  from  the 
solution  spectrograms  taken  during  primary  eclipse.  The 
first  actual  isolation  and  measurement  of  the  rotation 
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effect  in  stellar  spectra  was  done  by  Rossiter  (1924)  in 
3  Lyrae,  when  he  recognized  it  as  a  way  of  determining  the 
dimensions  of  an  eclipsing  system.  A  few  months  later 
McLaughlin  (1924)  substantiated  the  rotation  effect  in 
Algol,  and  subsequently  studied  it  in  a  Coronae  Borealis 
(McLaughlin  1934a),  and  in  X  Tauri,  6  Librae,  and  RZ 
Cassiopeiae  (McLaughlin  1934b) . 

Equations  for  the  calculation  of  the  rotation 
factor  for  uniformly  bright  spherical  stars  were  derived 
by  Petrie  (1938).  Equations  which  would  include  limb- 
darkening  were  presented  by  Kopal  (1942a;  1942b)  who,  in 
the  process,  introduced  and  developed  a  set  of  functions 
he  called  'associated  a-f unctions ' ,  which  are  generally 
applicable  to  an  analysis  of  other  effects  inherent  in 
observed  light  curves.  Up  to  this  point,  no  account  had 
been  taken  of  the  fact  that  the  axis  of  rotation  could  be 
inclined  from  the  axis  of  orbital  motion,  although  the 
above  authors  were  aware  that  it  would  have  an  effect. 

Then  Kopal  (1945)  made  a  fuller  analysis  which  took  this 
into  account,  and  subsequent  authors  (Hosokawa  1953;  Kita- 
mura  1953)  have  extended  and  refined  the  calculations. 


, 


CHAPTER  II 


OBSERVATIONAL  TECHNIQUES 

The  spectroscopic  eclipsing  binary  DE  Draconis  = 
HD  193964  has  been  observed  at  the  David  Dunlap  Observa¬ 
tory  (Toronto) ,  and  the  Dominion  Astrophysical  Observa¬ 
tory  (Victoria) ,  and  the  spectrograms  obtained  have 
yielded  the  radial  velocities  used  in  determining  the 
final  orbit.  This  chapter  discusses  the  methods  of  ob¬ 
serving  binary  systems,  in  particular  spectroscopic  and 
eclipsing  systems.  The  theory  behind  using  spectrograms 
to  measure  radial  velocities,  the  spectral  lines  that 
should  be  seen,  the  method  of  reduction,  and  some  sources 
of  error  are  also  discussed. 

Some  Basic  Theory 

The  basis  of  radial  velocity  investigations  is 
the  existence  of  the  well-known  Doppler-Fizeau  principle, 
which  is  applied  to  the  component  of  a  star's  motion 
along  the  line-of-sight ,  also  called  its  radial  velocity. 
The  change  in  wavelength  AA  is  proportional  to  the  true 
wavelength  A  and  the  ratio  of  the  relative  velocity  of 
the  star  and  the  observer,  V^_,  to  the  velocity  of  light, 
c: 

V  =  ^  (2.1) 

r  A 

The  velocity  is  positive  if  the  star  is  receding  from  the 
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observer  (the  observed  wavelength  is  longer) ,  and  is  neg¬ 
ative  if  the  star  is  approaching  the  observer  (shorter 
wavelength) .  Reliable  values  of  depend  on  the  accurate 
determination  of  the  true  wavelengths,  X,  and  the  precise 
measurement  of  AX. 

Most  stars  have  absorption  line  spectra,  where 
the  light  emitted  by  the  hotter  lower  photosphere  passes 
through  the  cooler  upper  photosphere.  The  wavelengths  of 
the  absorption  lines  indicate  which  elements  are  present 
in  the  atmosphere.  The  lines  from  the  star's  spectrum 
are  compared  with  those  of  a  known  source.  This  is  usu¬ 
ally  an  iron  arc  or  an  iron-argon  hollow  cathode  tube, 
whose  light  passes  over  the  same  path  in  the  spectrograph 
as  the  light  from  the  star.  This  comparison  spectrum  is 
recorded  on  each  side  of  the  stellar  spectrum. 

The  Spectrograph 

Most  radial  velocity  results  for  early-type  stars 
are  based  on  spectra  covering  wavelengths  from  3800  to 

O 

5000  A.  The  slit  of  the  spectrograph  is  illuminated  by 
the  image  of  the  star,  and  it  will  be  most  efficient  if 
the  slit  encompasses  the  complete  image.  However,  a  wide 
slit  reduces  the  spectral  purity,  which  is  measured  by 
the  width  of  the  projected  image  of  the  slit  at  the  re¬ 
cording  device  (photographic  plate) .  The  resolving  power 
of  the  spectrograph  is  usually  limited  by  the  slit  width 
and  not  by  the  dispersing  element. 


collimator  focussing 

element 


Fig.  2.  Schematic  diagram  of  a  spectrograph  (from 
Underhill  1966). 


The  effective  speed  of  a  spectrograph  is  dependent 
on  the  ratio  of  the  collimator's  diameter  to  the  diameter 
of  the  telescope  objective,  and  also  on  the  ratio  of  the 
collimator  focal  length  to  the  camera  focal  length.  The 
optimum  design  is  to  use  a  large  aperture,  long  focal 
length  collimator  and  a  short  focal  length  camera.  For 
high  resolution,  high  linear  dispersion,  gratings  are 
used  as  the  dispersing  element.  The  recording  device  is 
usually  a  photographic  plate  or  film,  with  the  requirement 
that  the  projected  slit  width  be  greater  than  the  size  of 
the  photographic  grain  (emulsion),  about  five  microns. 

The  plates  do  not  have  a  linear  response  to  various  light 
intensities,  so  an  intensity  calibration  spectrum  is  ex¬ 
posed  with  every  spectrogram  for  light  intensity  studies. 

A  comparison  with  the  reference  spectrum  external  to  the 
stellar  source  shows  the  displacement  of  the  lines  on  the 


stellar  spectrum  caused  by  the  Earth's  motion,  the  star's 
motion,  and  possibly  stellar  atmospheric  motion. 

Stellar  Spectrograms 

Radial  Velocity  Calculations 

Determining  the  wavelengths  from  spectrograms  in¬ 
volves  measuring  the  position  of  the  stellar  lines  with 
respect  to  the  comparison  spectrum  and  interpolating. 

This  procedure  involves  the  following  steps.  The  relation 
between  the  wavelength  and  the  positions  of  lines  in  the 
comparison  spectrum  is  called  the  dispersion  curve,  and 
is  given  by: 


X-X 

< 

s-s  =  — 
o  D 


(2.2) 


where  s  is  the  position  of  the  line  of  wavelength  X,  and 
X  ,  D,  and  sq  are  constants.  For  grating  spectrograms, 

D  is  the  linear  dispersion  in  angstroms  per  unit  setting. 

A  single  set  of  constants  for  a  given  spectrograph  can  be 

O 

used  over  a  range  of  about  200  A  by  making  use  of  a  cor¬ 
rection  curve.  The  settings  for  three  lines  of  known 
wavelengths,  one  at  each  end  and  one  in  the  middle  of  the 

reqion,  are  determined,  and  the  constants  X  ,  D,  and  s 
'  •  o  o 

are  calculated.  Then  the  positions  of  standard  wavelength 
in  the  comparison  and  typical  stellar  spectra  are  calcu¬ 
lated,  with  the  zero-point  of  the  scale  set  by  making  a 
certain  position  correspond  to  a  selected  wavelength  in 
the  comparison  spectrum;  for  this  study  the  zero-point 
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was  the  line  A4200  for  DAO  plates,  and  A3927  for  DDO 
plates . 

To  calculate  the  radial  velocity  of  the  star,  the 
following  procedure  was  used.  On  each  spectrogram,  the 
positions  of  both  the  standard  wavelengths  and  the  stellar 
spectral  lines  are  measured  after  setting  the  zero-point. 
At  DAO  the  comparison  spectrum  was  from  an  iron-argon 
tube  while  at  DDO  an  iron  arc  was  used.  The  measured 
positions  of  the  comparison  lines  will  not  be  the  same  as 
the  calculated  ones,  so  a  correction  curve  for  the  spec¬ 
trogram  is  drawn  of  the  difference  As  versus  wavelength. 

It  turns  out  to  be  approximately  parabolic  for  grating 
spectrographs.  Then  for  each  stellar  line,  the  correction 
to  the  measured  setting  at  that  wavelength  is  read  off 
this  curve  and  applied  to  the  measured  position.  The 
difference  between  the  measured  positions  of  the  stellar 
lines  and  the  calculated  positions  for  zero  velocity  are 
a  measure  of  the  star's  motion.  These  differences  are 
converted  to  a  velocity  by  the  factor  DVg,  where  D  is  as 
in  equation  2.2  (i.e.  the  slope  of  the  dispersion  curve 

at  that  wavelength)  ,  and  V  is  the  factor  -r-  : 

S  A 

V  =  DAs  y  =  AA  £  (2.3) 

r  a  A 

which  is  the  Doppler-Fizeau  principle.  Then,  for  each 
stellar  line  on  the  spectrogram,  a  is  obtained.  The 
average  of  these  is  taken  for  one  plate,  to  obtain  a 


calculated  for  that  phase  in  the  star’s  motion. 

However,  this  calculated  velocity  is  made  up  of 
the  star's  radial  velocity  with  reference  to  the  Sun, 
plus  the  component  of  the  Earth's  motion  in  the  line-of- 
sight.  The  latter  involves  three  components:  (1)  the 
Earth's  revolution  around  the  Sun;  (2)  the  Earth's  daily 
rotation  on  its  axis;  and  (3)  the  motion  of  the  Earth 
about  the  centre  of  mass  in  the  Earth-Moon  system.  The 
last  of  these  three  components  can  usually  be  omitted, 
being  of  the  order  of  ±0.014  km/s.  The  other  two  can  be 
computed  readily  (Underhill  1966),  with  the  correction 
for  the  Earth's  rotational  motion  varying  typically  from 
a  few  hundredths  to  a  few  tenths  of  a  km/s.  In  addition, 
if  a  prism  spectrograph  is  used,  a  correction  for  the 
curvature  of  the  spectral  lines  must  be  included  and  is 
characteristic  of  that  spectrograph.  With  most  grating 
spectrographs,  no  such  correction  is  required,  because 
the  images  of  a  long  slit  are  usually  straight  lines. 

In  order  to  monitor  the  stability  of  the  spectro¬ 
graph  and  of  his  measuring/reduction  techniques,  and  to 
permit  intercomparison  with  results  obtained  elsewhere, 
an  observer  normally  records  spectra  of  one  or  more 
'standard  velocity  stars'  each  night. 

Standard  Wavelengths 

It  is  most  important  that  accurate  zero-velocity 
wavelengths  be  determined  which  will  be  valid  for  a  given 
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dispersion,  applicable  over  a  range  of  spectral  types, 
and  representative  of  the  main  part  of  the  stellar  atmos¬ 
phere.  Most  stellar  spectrograms  of  high  dispersion  and 
high  resolution  show  many  lines  from  the  atoms  and  ions 
in  the  atmosphere.  In  early-type  spectra,  however,  there 
are  only  a  few  strong , relatively  broad  lines;  the  others 
are  weak  and  shallow.  Some  of  the  strong  lines  are  formed 
in  the  surface  layers  of  the  atmosphere,  and  their  posi¬ 
tion  in  the  spectrum  may  represent  the  peculiar  differen¬ 
tial  motions  of  these  layers,  rather  than  the  motion  of 
the  whole  star. 

Each  spectrograph,  with  its  own  optical  combina¬ 
tion  and  dispersion,  will  have  its  own  most  suitable  lines 
because  different  dispersions  cause  different  line-blend¬ 
ings.  The  procedure  of  calculating  appropriate  standard 
B-type  stellar  wavelengths  for  each  combination  of  spec¬ 
trograph-telescope-camera  involves  a  series  of  steps 
which  work  down  from  well-established  wavelengths  in  solar- 
type  spectra  to  well-behaved  B-type  control  stars.  First, 
a  set  of  solar  wavelengths  is  determined  whose  accuracies 
can  be  checked  by  comparing  the  radial  velocities  calcu¬ 
lated  using  these  values,  with  the  velocities  calculated 
from  positional  astronomy.  The  second  step  is  to  extend 
the  system  to  earlier  types  using  control  stars  of  types 
A  and  F  whose  true  radial  velocities  are  known.  These 
control  stars  either  belong  to  clusters  whose  motions  are 
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known  from  positional  calculations  and  radial  velocities 
of  solar-type  stars,  or  are  components  of  wide  binary 
systems,  having  a  companion  of  type  K5  to  F4,  and  small 
relative  orbital  motion.  The  third  step  is  to  extend  the 
set  of  wavelengths  to  B-type  stars  through  control  stars 
with  radial  velocities  known  by  other  means.  For  example, 
the  control  stars  are  members  of  visual  binaries  with  a 
companion  of  type  AO  or  later,  or  are  members  of  a  clus¬ 
ter  with  known  radial  velocity.  For  the  B-type  stars, 
standard  laboratory  wavelengths  are  initially  used,  and 
the  control  stars  indicate  the  suitability  of  each  line, 
its  stability  and  usefulness  for  main-sequence  or  super¬ 
giant  stars.  Difficulties  arise  from  the  diffuseness  of 
lines  in  B-type  spectra,  and  from  the  fact  that  one-half 
of  the  known  B-type  stars  are  spectroscopic  binaries  or 
have  large  scale  atmospheric  motions.  Most  of  the  suit¬ 
able  lines  are  due  to  Hel  and  hydrogen,  and  cover  the 
range  of  about  X3600  to  A4900. 

Interpolation  to  obtain  the  corrections  to  the 
measured  stellar  line  positions  also  calls  for  stable 
accurate  wavelengths  in  the  comparison  spectrum.  These 
are  called  secondary  standards,  and  Class  A  secondary 
standards  are  determined  to  an  accuracy  better  than 

O 

0.0001  A,  and  can  serve  as  substitutes  for  the  primary 

standard  which  is  a  line  from  Kr^  with  16057.80210^  in 

°  -10 

a  vacuum.  This  defines  1  A  as  10  meter  exactly. 


Class  B  secondary  standards  are  most  usually  used,  and 

O 

allow  the  measurement  of  lines  to  better  than  a  0.01  A 
accuracy.  The  wavelengths  have  an  accuracy  of  at  least 

O 

0.001  A.  They  can  be  lines  of  an  iron  arc  in  air  (which 
are  rather  broad)  or,  more  desirably,  iron  lines  emitted 
in  hollow  cathode  discharges  or  high-frequency  discharges 
using  Ne,  Ar,  or  He  as  the  carrier  gas.  The  standard 
wavelengths  used  in  reducing  the  DAO  plates,  obtained 
with  the  spectrograph-telescope  combination  described  in 
Chapter  3,  were  calculated  as  outlined  by  Batten  (1976). 


CHAPTER  III 


THE  SPECTROGRAPHIC  ORBIT  OF  DE  DRACONIS 

This  chapter  discusses  the  main  features  of  the 
spectroscopic  eclipsing  binary  system  DE  Draconis.  The 
methods  used  to  determine  the  period  and  final  orbital 
elements,  along  with  the  relevant  data,  are  outlined,  and 
the  results  are  compared  with  those  from  the  previous  re¬ 
search  efforts  of  Hube  (1973)  and  Willmarth  (1976).  In 
addition,  certain  elements  are  compared  with  those  ob¬ 
tained  from  photometric  observations.  The  photometric 
data  will  be  used  to  calculate  such  physical  parameters 
for  both  components  as  masses,  radii,  and  size  of  orbit. 
The  graph  of  the  final  velocity  curve  will  also  show  all 
observational  data  and  the  deviations  of  the  observed 
data  during  eclipse  from  the  calculated  values. 

Period  Determination 

In  calculating  the  elements  of  a  spectroscopic 
binary,  an  accurate  determination  of  the  period  is  the 
first  and  most  crucial  requisite.  This  can  be  very  dif¬ 
ficult  because  an  astronomer  is  restricted  by  the  weather, 
the  location  of  the  telescope,  and  the  allocation  of 
telescope  time  from  obtaining  many  observations  at  short, 
equal  intervals.  There  is  no  general  method  for  deter¬ 
mining  the  period,  and  all  methods  for  a  preliminary  value 
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are  based  on  trial  and  error.  With  computers,  the  number 
of  trials  made  is  increased  and  the  length  of  time  re¬ 
quired  to  complete  the  calculations  greatly  reduced.  An 
unknown  number  of  cycles  will  occur  between  different  ob¬ 
servations  at  recognizable  phases,  e.g.  maximum  or  minimum 
velocities,  so  that  short  periods  are  most  difficult  to 
recognize.  The  basic  method  is  to  pick  out  the  maxima 
and  minima  from  the  series  of  observations  at  different 
times.  Each  corresponding  pair  of  maxima  or  minima  will 
give  a  relation: 

t2  -  t  =  nP  (3.1) 

where  n  is  an  integer.  A  number  of  pairs  will  give  dif¬ 
ferent  multiples  of  the  period  P,  and  the  smallest  common 
to  all  is  taken  as  the  preliminary  value.  There  will 
usually  be  several  values,  multiples  of  P,  which  will 
satisfy  the  observations.  For  each  possible  value,  the 
phase  of  each  observation  is  calculated  and  the  observa¬ 
tions  are  plotted  to  determine  of  it  generates  an  accept¬ 
able  curve.  In  cases  where  it  is  difficult  to  choose, 
observations  must  be  made  on  certain  critical  dates  when 
the  alternative  periods  predict  very  different  velocities. 
However,  even  after  eliminating  all  the  identifiably  un¬ 
suitable  values,  the  remaining  period  may  not  be  the  true 
value.  One  serious  problem  is  that  in  the  case  of  spec¬ 
troscopic  binaries,  spurious  periods  may  arise  from  re¬ 
strictions  in  the  observing  program.  For  example,  with 
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any  one  telescope  and  location,  observations  are  often 
restricted  to  certain  periods  of  the  year.  Furthermore, 
these  observations  are  also  restricted  to  intervals  of 
about  one  sidereal  day  (±0.1  day),  because  measurements 
are  usually  taken  when  the  system  is  near  the  meridian. 
Tanner  (1948)  has  proposed  a  method  of  detecting  spurious 
periods  and  this  method  was  used  by  Hube  (1973)  in  his 
original  determination  of  the  period  of  DE  Draconis. 

Determination  of  Other  Orbital  Elements 

Once  the  preliminary  value  of  the  period  has  been 
found,  a  consideration  of  the  amplitude  and  shape  of  the 
best-fitting  curve  of  velocity  versus  phase  allows  the 
determination  of  the  elements  of  the  true  orbit  of  the 
star  with  respect  to  the  centre  of  gravity  of  the  system: 
e,  co ,  T,  K,  and  y.  In  addition,  the  projected  semi-major 
axis  asin(i)  and  the  mass  function  f (M)  can  be  calculated 
from  the  orbital  elements.  The  resulting  phase-velocity 
diagram  will  be  sinusoidal  unless  e  is  large.  The  ob¬ 
served, or  line-of-sight  component  of,  radial  velocity  with 
respect  to  the  Sun  consists  of  the  constant  velocity  of 
the  binary  system's  centre  of  mass,  y,  plus  the  variable 
orbital  velocity  of  the  star.  Applying  the  equations  of 
elliptic  motion  (Underhill  1966;  Aitken  1964),  the  final 
result  is: 

=  y  +  K[ecos(w)  +  cos  (v+oo)  ] 


(3.2) 
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where  v,  called  the  'true  anomaly' ,  is  the  angle  from 
periastron  to  the  star's  position  in  its  true  orbit  at 
any  time  t,  measured  in  the  direction  of  motion  (see 
figure  1).  Also,  the  element  K  can  be  expressed  in  terms 
of  P,  asin(i),  and  e. 

Preliminary  values  of  K,  e,  w,  v  (which  will  give 
T) ,  and  y  can  be  determined  through  direct  methods  (Under¬ 
hill  1966;  Aitken  1964)  using  the  initially  determined 
period  P  and  the  observed  velocity  data.  These  prelimin¬ 
ary  elements  are  then  refined  by  a  least-squares  solution, 
which  involves  the  simultaneous  adjustment  of  the  five 
parameters  (or  six  if  the  period  is  also  varied) ,  since 
vr  =  f(y*  K,  e,  oo,  T,  P)  .  For  each  time  of  observation 
the  difference  in  velocity  A (sometimes  denoted  by  O-C) 
is  the  observed  velocity  minus  the  calculated  velocity 
for  that  phase,  and  this  difference  is  produced  by  correc¬ 
tions  to  the  elements  of  Ay,  AK,  .... 

AVr  =  If  AT  +  M  AK  +  ••••  (3‘3) 

Assuming  all  the  changes,  or  corrections,  are  independent 
and  small  enough  to  include  only  first-order  terms,  the 
well-known  equation  of  condition  can  be  applied  to  each 
observation  (Lehmann-Filhds  1894;  Schlesinger  1908). 
Briefly  the  steps  are  the  following:  (1)  the  preliminary 
elements  are  used  to  calculate  the  radial  velocity  for 

8  f 

the  phase  of  each  observation;  (2)  the  co-efficients  t— , 
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...  of  the  differential  corrections  are  computed  for 
each  observation,  using  the  residuals  0-C  from  the  pre¬ 
liminary  values  in  step  (1);  (3)  normal  equations  are 

formed  and  solved  simultaneously  to  give  the  differential 
corrections  which  are  applied  to  the  initial  elements. 
These  steps  are  re-applied  any  number  of  times  to  obtain 
the  required  precision.  In  this  study,  the  values  were 
taken  as  the  final  orbital  elements  after  five  iterations, 
and  each  observation  was  given  equal  weight. 

Problems  arise  in  the  above  method  when  the  eccen¬ 
tricity  is  small.  For  a  circular  orbit,  e  =  0,  and  both 
T  and  to  have  no  significance.  When  e  is  small,  T  and  co 
are  almost  indeterminate  because  the  co-efficients  of  AT 
and  Ato  in  the  Lehmann-Filhes '  equations  approach  one,  and 
a  denominator  in  the  solution  of  the  normal  equations  is 
almost  zero.  Fixing  T  or  to  biases  the  solution  and  gives 
unrealistic  values  of  the  precision  of  the  elements. 

Sterne  (1941)  resolved  this  difficulty  by  using  the  date 
T^ ,  the  epoch  of  mean  longitude,  instead  of  T.  T^  is  the 
date  at  which  the  mean  longitude  (m  +  M)  is  zero,  where 
M  is  the  mean  anomaly,  defined  by: 

M  =  — -  (t  -  T)  (3.4) 

and  t  is  the  date  of  observation.  The  mean  longitude  at 
2  7T 

t  is  then  =—  (t-T  ) ,  and  T  remains  determinate  even  for 
P  o  o 

e  =  0.  When  e  is  small  but  non-zero,  a  modification  of 
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Schlesinger ' s  least-squares  method  to  correct  T  ,  instead 
of  T,  can  be  used.  But  if  e  is  nearly  zero  the  prelimin¬ 
ary  value  is  taken  to  be  e  =  0,  and  the  simple  modifica¬ 
tion  cannot  be  used.  Instead,  allowing  e  to  be  zero  ini¬ 
tially,  leads  to  new  equations  of  condition  which  are 
solved  by  least-squares  methods.  Thus,  modifications  to 
e  are  calculated  and  then  the  first  procedure  can  be 
applied.  The  second  method,  that  of  correcting  an  ini¬ 
tially  circular  orbit,  is  particularly  applicable  to  DE 
Draconis,  which  has  a  very  small  eccentricity. 

Previous  Orbital  Elements  of  DE  Draconis 

Although  DE  Draconis  =  HD  193964  =  HR  7792 

[a  (1950)  =  20h  18?8,  6  (1950)  =  +62°05',  m  =  +5.72,  B9V] 

v 

has  been  known  to  be  a  spectroscopic  binary  system  since 
the  early  part  of  this  century  (Plaskett  and  Young  1919), 
the  first  set  of  published  orbital  elements  was  presented 
by  Hube  (1973).  Hube  determined  a  preliminary  period 
using  28  spectrograms  from  the  David  Dunlap  Observatory 
and  four  velocities  from  Victoria  and,  because  the  radial 
velocity  curve  indicated  a  low  eccentricity,  initially 
determined  the  best-fitting  circular  orbit.  Employing 
Sterne's  (1941)  method,  corrections  to  the  circular  orbit¬ 
al  elements  were  calculated,  and  the  corrected  values 
were  then  used  as  the  preliminary  quantities  in  the  appli¬ 
cation  of  Schlesinger ' s  method  (1908).  Hube's  final 
elements  are  given  in  table  1.  Hube  used  Tanner's  method 
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(1948)  to  detect  any  spurious  periods  but  obtained  nega¬ 
tive  results.  Willmarth  (1976)  performed  an  independent 
spectroscopic  investigation  of  DE  Draconis  using  eleven 
velocities  from  earlier  investigations  plus  fifteen  new 
velocities  from  the  Ritter  Observatory.  The  Willmarth 
results  are  also  shown  in  table  1. 

Comparison  of  the  two  sets  of  elements  shows 
large  differences  in  the  systemic  velocity  and  the  semi¬ 
amplitude.  Hube  (1976)  obtained  another  solution  using 
all  the  published  velocities  to  that  date,  in  which  the 
period  was  5^298095  ±  0.000025.  Hube  pointed  out,  as  he 
had  in  his  original  paper  (Hube  1973),  the  need  to  clarify 
the  scatter  of  the  observed  radial  velocities  near  the 
time  of  conjunction,  which  was  larger  than  would  have 
been  expected  from  internal  errors. 

Fuertig  (1975)  found  DE  Draconis  to  be  an  eclip¬ 
sing  binary  and,  using  Hube ' s  (1973)  period  and  working 
with  Mueninger  (Fuertig  and  Mueninger  1976),  made  further 
observations  of  the  primary  minimum  (a  transit) .  The 
preliminary  photometric  elements  derived  by  Fuertig  and 
Mueninger  are  presented  in  table  2.  Fuertig  and  Mueninger 
concluded  from  these  elements  that  the  invisible  companion 
is  a  solar-type  star,  as  expected  from  Hube ' s  (1973)  spec¬ 
troscopic  observations.  No  secondary  minimum  was  observed. 


THE  ORBITAL  ELEMENTS 
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TABLE  2 

PHOTOMETRIC  ELEMENTS 


Co-efficient  of 

limb- 

-darkening 

of  large  star 

X 

g 

0.5 

Ratio  of  radii 

k 

0.35 

Radius  of  large 

star 

(a  =  1.0) 

r 

g 

0.157 

Radius  of  small 

star 

(a  =  1.0) 

r 

s 

0.055 

Inclination  of  < 

orbit 

plane 

i  (°) 

86.85 

Eccentricity 

e 

0 

Period 

P  (days) 

5.2984 

New  Observations 

During  1976  and  1977  a  total  of  87  useful  new 
spectrograms  of  DE  Draconis  were  obtained:  one  at  the 
David  Dunlap  Observatory  (DDO)  and  the  remainder  at  the 
Dominion  Astrophysical  Observatory  (DAO) .  These  spectro¬ 
grams  were  augmented  with  18  higher  dispersion  plates 
from  the  DDO  that  Hube  gathered  between  1970  and  1972 
for  his  original  calculations.  The  DDO  observations  were 
obtained  with  the  grating  spectrograph  on  the  1.88m  re- 

O 

flector  telescope,  giving  a  dispersion  of  12  A/mm.  Typi¬ 
cally,  on  these  plates  between  six  and  twelve  lines  could 
be  measured,  all  plates  having  similar  internal  probable 
errors.  The  DAO  plates  were  taken  with  the  21121  spectro¬ 
graph  at  the  Cassegrain  focus  of  the  1.8m  telescope, 


ji 


38 


giving  a  dispersion  of  15  A/mm.  On  these  spectrograms 
usually  from  six  to  ten  lines  could  be  measured,  and  all 
plates  had  internal  errors  similar  to  those  from  the  DDO 
plates.  Thus,  all  velocities  were  given  equal  weight. 
The  DAO  plates  were  measured  with  the  oscilloscope  com¬ 
parator  at  Victoria  and  reduced  using  the  standard  wave¬ 
lengths  for  this  system  as  published  by  Batten  (1976). 
These  observations  are  documented  in  Appendix  A. 

Standard  velocity  stars  were  observed  on  each 
night  as  well,  and  the  calculated  velocities  of  these 
stars  indicated  that  a  correction  of  +1.7  km/s  should  be 
added  to  the  velocities  obtained  from  the  DAO  plates. 

The  observations  in  Appendix  A  from  the  DAO  include  this 
correction  factor. 


New  Orbital  Elements 

In  deriving  orbital  elements,  all  those  velocities 

cl  cl 

obtained  during  eclipse  (phase  about  1.7  to  2.0)  were 
omitted,  resulting  in  52  of  the  total  105  observations 
being  used.  The  method  of  Lehmann-Filhds  was  applied 
over  five  iterations  to  obtain  differential  corrections 
to  six  elements.  The  preliminary  values  used  were  those 
of  Hube  (1973) .  The  resultant  period  was  then  fixed,  and 
the  remaining  five  elements  were  calculated  from  five 
iterations  of  the  Lehmann-Filhds  equations  of  condition 
for  five  elements.  The  final  values,  the  mean  errors, 
and  Hube’s  (1973)  original  solution  are  listed  in  table  3. 


FINAL  ORBITAL  ELEMENTS 
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In  figures  3  and  4,  the  observations — the  52  out¬ 
side  eclipse  and  the  53  during  eclipse — have  been  plotted 
along  with  the  final  velocity  curve  from  the  Lehmann- 
Filhds  solution.  It  shows  clearly  the  rotation  effect 
during  the  time  of  eclipse  when  the  smaller  (invisible) 
secondary  transits  in  front  of  the  larger,  brighter  prim¬ 
ary  component. 

The  apparent  changes  in  the  orbital  elements 
(see  table  3)  are  probably  not  real;  the  earlier  velocity 
measurements  were  not  as  carefully  controlled  with  obser¬ 
vations  of  standards  as  were  the  later  observations. 
Comparison  of  the  two  sets  of  elements  shows  that  the 
calculated  systemic  velocity  has  apparently  decreased,  as 
has  the  eccentricity.  The  longitude  of  periastron  has 
increased  but  because  of  the  nearly  circular  orbit,  this 
quantity  will  be  very  uncertain.  The  photometric  data  of 
Fuertig  and  Mueninger  (1976)  yielded  a  period  of  5^2984 

and  a  duration  of  the  eclipse  of  about  0.06  of  a  period, 

d  d 

or  0.32.  The  times  of  conjunction  should  occur  0.16 

after  the  start  of  the  eclipse.  From  the  spectroscopic 

observations  in  figure  4,  the  eclipse  appears  to  last 

from  roughly  phase  1^68  to  phase  1^98,  or  about  0^30. 

The  time  of  conjunction,  when  the  curve  goes  through  the 

systemic  velocity  line,  comes  at  about  phase  1.83,  or 

d  d 

0715  after  the  start  of  the  eclipse  at  phase  1.68.  Thus, 

the  photometric  results  and  the  spectroscopic  calculations 


Fig.  3.  Final  velocity  curve  based  on  the  52  plates 
obtained  outside  eclipse.  (See  Appendix  A  for  data.) 
Crosses  represent  velocities  from  DDO  (Toronto) ;  open 
circles  indicate  DAO  (Victoria)  values. 
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Fig.  4.  Observed  velocities  within  eclipse.  The  53  points  fit  in  the  blank  section 
of  fig.  3.  Crosses  represent  DDO  velocities;  open  circles,  DAO  plates. 
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show  good  agreement. 

No  sure  double-lined  spectra  have  thus  far  been 
obtained,  but  there  are  indications  on  some  plates  of  line¬ 
splitting.  Currently,  spectrograms  of  DE  Draconis  are 
being  obtained  during  primary  eclipse  times  with  a  higher 
dispersion  spectrograph  on  the  48-inch  telescope  at  DAO. 
Outside  of  eclipse,  the  brighter  primary  overwhelms  the 
secondary,  but  during  the  start  and  end  of  the  eclipse, 
when  part  of  the  primary  is  blocked  off,  it  is  hoped  that 
a  higher  dispersion  will  give  plates  with  distinct  line 
structure  changes. 

Photometric  Data  and  Physical  Parameters 

Crawford,  Barnes,  Golson,  and  Hube  (1973)  have 
published  results  of  uvbyB  photometry  for  DE  Draconis  and 
obtained  the  following  parameters:  3  =  2.824;  b-y  =  -0.017; 
m^  =  0.123;  c^  =  0.873;  V  =  +5.72.  These  colour  para¬ 
meters  have  been  used  to  calculate  some  further  physical 
characteristics  for  DE  Draconis.  The  details  may  be  found 
in  Appendix  B,  but  a  brief  description  of  the  methods  and 
results  is  presented  here.  6  is  the  strength  of  the  H8  line. 

It  should  be  noted  that  the  above  data  refer  to 
the  integrated  light  from  both  components.  Golay  (1972) 
pointed  out  that  binarity  has  an  effect  which  is  photome¬ 
trically  similar  to  a  small  change  in  spectral  type. 
Theoretically  (McNamara  1966;  Golay  1974)  the  effect  of 
adding  companions  of  various  masses  to  primary  stars  of 
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different  spectral  types  displaces  the  colour-luminosity 
diagram  (M  versus  B-V)  above  and  approximately  parallel 
to  the  main-sequence  line  appropriate  for  single  stars  of 
the  primary  spectral  types.  In  a  colour-colour  diagram 
(U-B  versus  B-V) ,  the  linear  part  (primary  types  BO  to 
about  A5)  remains  very  close  to  the  single  star  sequence 
i.e.  with  the  inclusion  of  a  companion,  the  B-V  value  re¬ 
mains  essentially  the  same,  while  the  U-B  value  increases. 
In  the  curved  part  (primary  types  FO  to  G5) ,  the  change 
in  the  B-V  value  is  more  noticeable,  while  the  change  in 
U-B  is  less  so.  For  DE  Draconis,  with  a  primary  of  type 
B9V  and  mass  ratio  about  3.7  (see  Appendix  B) ,  the  theo¬ 
retical  effect  of  the  secondary  on  the  absolute  magnitude 
and  colours  assigned  to  the  primary  would  seem  to  be  un- 
noticeable,  judging  from  the  diagrams  and  tables  given 
by  Golay  (1974,  pp.  91-93).  Also,  the  following  rough 
calculations  will  show  that  in  the  system  the  light  from 
the  secondary  contributes  a  negligible  amount  to  the 
total,  and  therefore  the  results  from  the  uvby3  photometry 
can  be  assumed  to  give  accurate  information  on  the  primary 
alone . 

Assuming  the  two  components  of  DE  Draconis  to  be 
B9V  and  G2V,  and  using  colour  and  magnitude  values  appro¬ 
priate  on  average  for  these  spectral  types,  the  fraction 
of  light  contributed  by  the  secondary  at  different  wave¬ 
lengths  can  be  calculated  from  the  equation: 


' 


46 


and 


=  -2.5  log(L^) 

log  AA  =  0.4  (Mu  -  M2a) 
L1X 


(3.5) 

(3.6) 


where  X  is  any  wavelength,  and  subscripts  1  and  2  refer 
to  the  primary  and  secondary  respectively.  The  results 
for  UBV  wavelengths  are  shown  in  table  4 . 


TABLE  4 

FRACTIONAL  LIGHT  FROM  THE  SECONDARY 
(PRIMARY  =  1.0) 


A 

M 

B 

M 

U 

L2V 

L2B 

L2U 

G2V 

+  4.7 

+  5.34 

+  5.50 

0.025 

0.013  0.009 

B9V 

+  0.7 

+  0.64 

+  0.45 

From 

the  above 

results , 

it  is 

obvious  that 

the  contribu- 

tion 

from  the 

secondary 

at  any 

wavelength 

in  the  visible 

region  can  be  safely  ignored,  and  the  uvby3  data  can  be 
assumed  to  apply  to  the  primary  star  alone  with  little  or 
no  loss  of  accuracy. 

It  is  interesting  to  note  here  the  effect  that 
the  axial  rotation  of  the  B9V  component  could  have  on  the 
colour  parameters.  Kraft  (1970)  investigated  the  in¬ 
fluence  of  rotation  on  luminosities  and  temperatures  of 
main-sequence  stars.  A  rotating  star  which  is  distorted 
from  a  spherical  shape  has  a  smaller  surface  gravity  at 
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the  equator  than  at  the  poles,  so  the  surface  brightness 
and  local  effective  temperature  are  lower  at  the  equator 
than  at  the  poles.  Therefore,  mean  luminosity  and  temp¬ 
erature  are  a  function  of  the  angle  between  the  line-of- 
sight  and  the  axis  of  rotation.  Comparing  a  rotator  and 
a  non-rotator  of  the  same  mass  should  in  principle  show 
that  the  rotator  has  a  slightly  different  luminosity, 
since  the  gravitational  potential  and  therefore  the  cen¬ 
tral  pressure  and  temperature  are  changed. 

Theoretical  results  (Kraft  1970;  Roxburgh  and 
Strittmatter  1965)  indicate  that  a  rotator  is  cooler  than 
a  non-rotator,  regardless  of  aspect,  and  also  that  a 
rotator  is  brighter  (if  seen  pole-on)  or  fainter  (if  seen 
equator-on)  than  a  non-rotator  of  the  same  mass.  How¬ 
ever,  these  effects  are  very  small,  amounting  to  at  most 
a  few  hundred  degrees  in  T  and  a  few  tenths  of  a  mag¬ 
nitude  in  .  The  main  sequence  for  rotating  stars  is 

expected  to  lie  to  the  right  of  the  non-rotator  main  se- 

2 

quence,  with  the  displacement  depending  mainly  on  V  (the 
rotational  velocity  squared) ,  and  not  so  much  on  i  (the 
angle  between  the  axis  of  rotation  and  the  line-of-sight ) . 
(See  Strittmatter  and  Sargent  1966;  Hardorp  and  Stritt¬ 
matter  1968).  For  early-type  stars  in  the  U-B  versus  B-V 
diagram,  rotation  projects  stars  along  the  main  sequence 
for  non-rotators,  rather  than  at  an  angle  to  it.  A  simi¬ 
lar  effect  occurs  in  the  c^  versus  (b-y)  diagram,  indi- 
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eating  that  rotation  does  not  apparently  effect  the  colour- 
colour  diagrams  for  early-type  stars. 

Using  a  recent  calibration  of  the  uvby$  system 
(Crawford  1975)  ,  the  value  of  c^  was  corrected  for  red¬ 
dening  to  obtain  the  unreddened  parameter  cq  for  DE  Dra- 

conis.  From  Crawford’s  (1978)  calibration  of  c  versus 

o 

3  (Zero  Age  Main  Sequence  =  ZAMS) ,  the  value  of  3  for  a 
ZAMS  star  was  determined.  Crawford  (1978)  also  calibrated 
My  (ZAMS)  versus  3  (ZAMS)  ,  and  gave  a  correction  term  for 
stars  lying  above  the  ZAMS  (evolved  stars) .  A  value  for 
My (ZAMS)  was  obtained  for  DE  Draconis  and  corrected  for 
evolutionary  effects,  since  the  3  (observed)  differed  from 
the  3  (ZAMS).  The  absolute  visual  magnitude  for  DE  Dra¬ 
conis  was  calculated  to  be  My  =  +0.75,  which  compares 
well  with  the  MK  value  of  +0.6,  and  also  with  the  value 
adopted  by  McNamara  (1966)  of  +0.7. 

Crawford's  (1975)  description  of  the  uvby3  system 
also  allows  the  apparent  visual  magnitude  to  be  corrected 
for  reddening,  and  from  the  distance  modulus  (rriy  -  My), 
the  distance  to  this  system  was  found  to  be  about  96 
parsecs . 

Using  Crawford's  (1978)  relation  between  cq  and 

(U-B)  ,  calculations  for  0  were  made.  Several 

°  etr  eff 

different  formulae  (Heintze  1973)  yielded  similar  values 
for  Teffr  and  the  finally  adopted  average  value  of 
11 , 300°K  agrees  well  with  values  for  spectral  type  B9v 
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from  other  effective  temperature  versus  spectral  type 
scales  (e.g.  McNamara  1966;  Hack  and  Struve  1969)  . 

From  the  formula  given  in  Allen  (1973) ,  the  bolo- 
metric  correction  was  calculated,  which  gave  an  absolute 
bolometric  magnitude  =  +0.19.  The  radius  was  calcu¬ 

lated  from  two  different  formulae  (Allen  1973;  Kopal  1959) 
and  the  average  radius  of  DE  Draconis  is  2.2  solar  radii. 
The  mass  was  also  calculated  using  various  mass-luminosity 
relationships  (Batten  1968;  Allen  1973;  Petrie  1950; 
Harris,  Strand,  and  Worley  1963;  Kopal  1959),  and  the 
average  mass  is  2.7  solar  masses.  These  values  of  mass 
and  radius  agree  well  with  those  expected  for  a  typical 
B9V  star  (see  Underhill  1966) . 

Further  results  can  be  found  for  the  secondary 
component  using  the  mass  function,  projected  semi-major 
axis,  and  the  photometric  elements  of  Fuertig  and  Muenin- 
ger  (1976).  From  the  inclination  i  =  86°85,  the  calcu¬ 
lated  mass  of  the  primary,  and  the  mass  function,  the 
secondary  is  found  to  be  about  one  solar  mass.  This  is 
consistent  with  the  expectations  from  the  photometric 
data  of  Fuertig  and  Mueninger  (1976)  for  a  solar-type 
secondary,  since  no  secondary  minimum  had  been  observed. 

It  is  also  in  agreement  with  the  spectroscopic  observa¬ 
tions,  where  no  lines  of  the  secondary  were  visible.  The 
radius  of  the  relative  orbit,  a  =  a^  +  a2 r  was  also  cal¬ 
culated,  and  from  the  photometrically  determined  frac- 
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tional  radii,  the  absolute  radii  were  found  to  be  3.1  and 
1.1  solar  radii  for  the  primary  and  secondary,  respective 


The  values  of  the  radii  above  differ  from  those 
obtained  using  the  calculated  bolometric  absolute  magni¬ 
tude,  but  still  conform  within  the  limits  expected  for 
these  spectral  types.  Also,  the  calculation  of  ^  is 
liable  to  be  very  approximate  because  there  appear  to  be 
many  differences  of  opinion  regarding  the  correct  effec¬ 
tive  temperature  calculation,  and  mass-luminosity  rela¬ 
tionship.  The  parameters  derived  above  and  some  compari¬ 
son  data  are  summarised  in  table  5.  The  values  of  ^ 
shown  in  brackets  were  computed  using  the  calculated 

The  adopted  values  of  3.1  and  1.1  solar  radii  were 
determined  using  photometric  data,  and  the  values  for 

a  9  were  calculated  using  these  adopted  radii. 

1 ,  ^ 


. 


PHYSICAL  PARAMETERS  OF  DE  DRACONIS 
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CHAPTER  IV 


THE  ROTATION  EFFECT 

In  this  chapter  we  will  discuss  the  distortion  of 
a  radial  velocity  curve  of  a  binary  system  due  to  the  eclipse 
of  an  axially  rotating  star.  Some  straightforward  effects 
which  can  be  expected  and  the  general  theory  required  to  cal¬ 
culate  these  effects  under  certain  assumptions  will  be  pre¬ 
sented.  The  results  of  calculations  done  for  DE  Dra,  and  the 
implications  for  certain  physical  dimensions  of  this  system 
will  be  presented. 


General  Information 

Probably  all  stars  rotate,  although  to  varying  de¬ 
grees.  For  the  more  rapidly  rotating  stars  the  effect  shows 
up  in  the  spectral  lines,  and  the  rotational  velocity  Vesina 
can  be  calculated  if  Vasina  is  greater  than  about  25  km/s. 
(Here  a  is  the  angle  between  the  axis  of  rotation  and  the 
line-of-sight ) .  The  total  light  received  from  a  star  is 
made  up  of  contributions  from  the  visible  part  of  the  hemi¬ 
sphere.  For  a  single  rotating  star  one  limb  is  approaching 
the  observer,  and  the  other  limb  is  receding.  The  light  from 
the  approaching  limb  thus  shows  absorption  lines  displaced 
toward  the  violet,  while  for  the  other  limb  they  are  dis- 
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placed  toward  the  red.  The  overall  effect  is  a  set  of  broad¬ 
ened  absorption  lines,  the  widths  depending  on  the  star's 
rate  of  rotation  and  orientation.  For  a  rapid  rotator  the 
lines  are  also  shallower  Cless  intense)  than  for  a  slower  ro¬ 
tator  because  the  absorption  is  spread  over  a  greater  region 
of  the  spectrum.  The  centres  of  the  lines  are  not  displaced 
by  the  rotation  when  the  whole  stellar  disc  is  observed.  If 
the  rotating  star  happens  to  be  a  component  of  an  eclipsing 
binary  system,  its  spectral  lines  are  asymmetrically  broad¬ 
ened  during  the  eclipse.  The  observed  velocity  is  not  only 
the  orbital  velocity  but  also  includes  a  contribution  from 
the  rotation.  If  (as  is  almost  always  the  case)  the  direc¬ 
tions  of  axial  rotation  and  orbital  revolution  are  the  same, 
then  at  the  start  of  eclipse  the  other  component  will  obscure 
that  limb  of  the  rotating  star  which  is  approaching  the  ob¬ 
server.  The  receding  limb  (shift  to  the  red)  contributes 
and  the  apparent  velocity  is  more  positive  than  the  calcu¬ 
lated  orbital  velocity.  At  the  end  of  eclipse  only  the  ap¬ 
proaching  limb  is  visible  (shift  to  the  violet)  and  observed 
velocities  are  more  negative.  If  the  eclipse  is  total,  then 
the  star  is  not  visible  for  the  duration  of  totality,  but 
just  prior  to  and  after  totality  a  thin  crescent  of  the  ro¬ 
tating  star  is  seen.  In  this  case  even  a  modest  rotational 
velocity  can  produce  large  departures  from  the  expected 
values.  For  an  annular  eclipse  (as  in  DE  Dra)  a  ring  of  the 
eclipsed  rotating  star  is  always  visible  and  light  from  one 
side  is  at  least  partly  balanced  by  light  from  the  other. 
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The  effect  is  thus  reduced,  however  a  high  rotational  velo¬ 
city  will  produce  a  small  effect.  At  the  middle  of  an  annu¬ 
lar  eclipse  the  opposite  sides  of  the  ring  exactly  balance 
and  the  net  rotational  effect  is  zero.  The  central  portion 
of  the  star,  which  has  no  Doppler  effect  due  to  rotation  but 
produces  a  deeper  centre  of  the  absorption  line,  is  eclipsed. 
At  this  point  the  absorption  lines  are  shallow,  and  some  may 
have  a  central  peak  as  if  they  were  doubled  lines.  Figure  5 
illustrates  the  eclipse  configurations  for  a  transit  and  the 
effect  on  the  radial  velocity  curve.  This  situation  is 
typified  by  DE  Dra. 


Fig.  5  . 


The  configuration  for  a  transit 
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The  observed  rotational  velocities  are  the  projec¬ 
tions  on  the  line-of-sight  of  the  true  equatorial  velocities: 

Vobs=  vesina*  If  the  assumption  is  made  that  the  axes  of 

_  4 _ 

rotation  are  randomly  oriented  in  space  then:  ’ Vasina. 

This  assumption  is  justified  by  agreement  of  the  calculated 
distribution  of  Vasina  assuming  random  orientation  with  the 
observed  distribution  of  Chandrasekhar  and  Munch  (1950) . 

Calculations  of  the  rotational  velocities  from  the 
spectral  line  shapes  have  been  made  for  stars  belonging  to 
all  spectral  types  and  luminosities.  The  results  show  that 
the  rotational  velocities  depend  on  the  spectral  type  and 
luminosity.  Hack  and  Struve  (1969)  present  the  results  of 
many  investigations.  In  particular  it  is  found  that  the 
early-type  stars  have  the  largest  observed  rotational  velo¬ 
cities.  For  main-sequence  stars,  the  values  are  greatest 
between  types  05  and  B5:  -  160  km/s  (Hack  and  Struve  1969, 

fig.  4b).  For  B9V  stars  -  130  km/s.  However  there  are 

stars,  notably  Be  stars  and  shell  stars,  with  -  340  km/s 

and  520  km/s,  respectively.  Slettebak  (1966)  calculated  the 
break-up  velocities — when  the  centrifugal  acceleration  be¬ 
comes  equal  to  the  gravitational  acceleration--based  on  cer¬ 
tain  main-sequence  rotating  star  models.  He  found  that  these 
break-up  velocities  are  higher  than  the  greatest  observed 
rotational  velocities:  for  a  B9V  star,  an  estimate  from  his 

curve  shows  V  (break-up)  -  450±25  km/s  while  V  (max.  obs.)  - 

0  ^ 

350±25  km/s.  The  difference  between  the  two  curves  is  a  min¬ 
imum  for  the  earliest  B  stars  and  increases  with  advancing 


.  -I 
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spectral  type. 

The  rotation  of  stars  which  are  members  of  binary 
systems  has  also  been  investigated.  Plavec  (1970)  points 
out  that  Plaut  (1959)  concluded  that  for  systems  with  orbi¬ 
tal  periods  less  than  10  days,  the  periods  of  rotation  and 
revolution  are  usually  the  same,  while  for  longer  periods 
there  is  no  synchronism  and  the  stars  rotate  like  single 
objects.  Plavec  (1970)  observed  a  definite  tendency  for 
components  of  binary  systems  to  rotate  more  slowly  than  the 
average  for  single  stars  of  the  same  type.  For  example,  an 
estimate  from  his  plot  yields  [single)  ^  190  km/s  but 
Vg  (binary)  -  80  km/s.  He  also  plotted  actual  rotational 

velocities  versus  calculated  synchronous  velocities  for  close 

2 tt  (radius) 


binary  systems,  V 


,  and  found  that  synchronism 


syn  period 

applies  at  least  approximately  in  most  cases.  There  were  no 
cases  where  the  component  rotated  much  more  slowly  than  was 
required  for  synchronism,  but  a  few  stars  rotated  much  more 
rapidly.  Plavec  also  noted  that  in  spectroscopic  binaries 
of  short  or  moderate  period,  the  secondaries,  if  less  than 
1M  ,  can  always  be  expected  to  rotate  in  synchronism  with 
the  orbital  motion. 

The  observation  that  close  binaries  generally  ro¬ 
tate  more  slowly  than  single  stars  is  supported  by  Jaschek 
(1970).  He  found  for  B  stars  that  for  members  of  spectro¬ 
scopic  binaries  the  largest  percentage  has  Vasina  in  the 
range  0  to  100  km/s,  while  for  normal  single  stars  the 
range  with  the  highest  percentage  is  from  100  to  200  km/s. 
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As  will  be  seen  later,  the  rotational  velocity  calculated  for 
DE  Dra  lies  at  the  upper  end  of  the  range  expected  for  single 
stars,  and  is  also  much  greater  than  the  synchronous  velocity. 
This  contradicts  the  above  general  observations.  However 
Nariai  (1971)  investigated  various  published  rotational  vel¬ 
ocities  for  double-lined  spectroscopic  and  eclipsing  binaries 
and  found  that  for  normal  main-sequence  stars  as  components 
of  double-  and  single-line  binaries  the  deviation  from  syn¬ 
chronism  is  large. 


Preliminary  Theory 

The  displacement  of  a  spectral  line  from  its  zero 
velocity  position  involves  contributions  to  the  line-of- 
sight  component  from  the  visible  area  of  the  star.  The 
line-of-sight  components  of  velocity  of  points  on  a  spheri¬ 
cal  star  rotating  as  a  rigid  body  are  the  same  for  all  points 
on  a  line  on  the  apparent  (projected)  disc  parallel  to  the 
projected  axis  of  rotation.  Assuming  the  axis  of  rotation 
is  perpendicular  to  the  orbital  plane  (a  =  i) ,  the  following 


Let  dA  =  area  of  a  strip  of  the 

disc  parallel  to  the  pro¬ 
jected  axis  of  rotation 
a  =  distance  of  the  strip 
from  the  axis 

i  =  inclination  of  the  orb¬ 
ital  plane  to  the  plane 
of  the  sky 
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Ve  =  rotational  velocity  of  a  point  on  the  equator 

Vr  =  line-of-sight  component  of  the  velocity  of  strip  dA 

Then  V  =  V  sini  (4.1) 

re 

and  the  deviation  of  the  radial  velocity  from  purely  ellip¬ 
tic  motion  of  the  visible  area  of  the  star  is: 


V 


obs 


=  V  sini 
e 


/  adA 
/dA 


(4.2) 


where  the  integration  is  over  the  visible  portion  only, 

that  the  'rotation  factor'  F  is: 

p  =  /adA 
r  “  /dA 


so 


(4.3) 


Thus  F  is  the  mean  value  of  the  rotational  component  in  the 
line-of-sight  averaged  over  the  visible  portion  of  the  star. 


Petrie's  Analysis 

Petrie  (1938)  performed  the  above  integration  for 
an  eclipse  of  a  spherical,  uniformly  bright,  rigidly  rotating 
star  in  a  circular  orbit.  The  origin  of  the  co-ordinate 
system  is  the  centre  of  the  brighter  star  and  the  X-axis  is 
along  the  line  joining  the  centres  of  the  projected  discs. 
Petrie's  basic  equation  is: 


//| dxdy 
// dxdy 


(4.4) 


where  s  =  distance  of  a  point,  on  the  projected  disc,  from 
the  projected  axis  of  rotation 
and  k  =  radius  of  the  small  star. 

Petrie  evaluated  the  equation  for  three  cases:  (1)  a  small 
bright  star  eclipsed  by  a  larger  fainter  star;  (2)  a  large 
bright  star  eclipsed  by  its  smaller  fainter  companion;  and 


. 
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(3)  eclipses  of  ellipsoidal  stars.  DE  Dra  is  an  example  of 
case  (2) . 

If  the  radius  of  the  larger  brighter  star  (whose 
centre  is  the  origin)  is  unity,  the  projected  distance  be¬ 
tween  the  centres  is  given  by: 

1  / - 2 - ‘ — 2 - T~ 

6  =  —  /cos  i  +  sin  isin  0  (4.5) 

r  £ 

where  r^  =  radius  of  the  large  star  with  respect  to  the  rad¬ 
ius  of  the  relative  orbit 
0  =  longitude  from  mid-eclipse 


To  observer 

Fig.  7 .  Projected  configuration  of  eclipse  (from  Petrie 
1938) 


' 


The  distance  s  of  a  point  on  the  primary  from  the  projected 
axis  is  then: 
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s  =  xcos4>  -  ysin<f> 


(4.6) 


where 


COSCj)  = 


sinG 


6r 


and 


sincf)  = 


cosicosS 


i 


6  r 


(4.7) 


i 


The  boundaries  of  the  projected  discs  are  circles  defined  by: 
2  2 

x  +  y  =1  larger  brighter  star  (4.8) 

2  2  2 

(x-6)  +  y  =  k  smaller  fainter  star  (4.9) 


The  points  of  intersection  of  the  two  circles  have  co-ordin¬ 
ates  : 


x 


1 


y 


l 


1+6' 


26 


(1  +  6' 


2  2 
-  kV 


46' 


(4.10) 

(4.11) 


The  equation  for  F  becomes: 


_  rr  //(xcosfj)  -  ysinc}))dxdy 


F  =  k 


(4.12) 


f  fdxdy 


where  in  the  numerator  the  x-integration  for  the  visible  por¬ 
tion  is  done  first  from  the  limb  of  the  brighter  to  the 
eclipsing  limb  of  the  fainter.  The  y-integration  then  sums 
these  strips  from  one  intersection  point  to  the  other.  The 
parts  of  the  disc  above  and  below  this  cancel  each  other, 
and  thus  appear  only  in  the  denominator  which  must  sum  over 
the  whole  visible  part  of  the  star.  Therefore  three  cases 
arise  from  this  method  of  performing  the  integration  as  il¬ 
lustrated  in  the  following  Figures  8a,  8b,  and  8c. 


■ 


' 
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Fig.  8.  Configurations  for  Petrie's  three  cases 


The  corresponding  formulae  are: 
Case  1  (Fig.  8a):  <_  6 


F  = 


2  2 

f.  ,  r  ,1-6  -k  v  , 

-<5coscf)  i  ( - ^ - )y1  +  y1 


A2-y  2  +  k2sin 


7T  +  2y16  -  yx 


yi  -  yi 


-1  ,.2  .  -1 


X:^-y?  -  sin  J'y1  -  k~sin  ‘‘"/l-, 

1  1  (F"} 


(4.13) 


Case  2  (Fig.  8b) :  1  >  x^  >  6 


2  ,  2 


*  ,  r  ,1-6 “-k^ x  ^  ,2  J71  2  2  .  -l.yl., 

-5coscf>  l( - ^ - )y^  +  Trk  -  y^/k  -y  -  k  sm  (— )  } 


F  = 


IT  (1-k2)  +  2y16-y1/l-y^  +y. 


/7l  2  .  -1  ^  ,2  .  -1/1. 

/k  -y  -  sm  y^  +  k  sm  {—) 


(4.14) 


Case  3  (Fig.  8c) :  6  <  1-k 


_  _  -k  6coscJ) 

F  2 
1-k 


(4.15) 


The  results  of  applying  Petrie's  analysis  to  the 
observational  points  from  eclipse  data  of  DE  Dra  for 

several  different  values  of  the  rotational  velocity  are  shown 
in  Table  6  and  Figure  9 .  There  is  a  large  amount  of  scatter 


' 


ROTATIONAL  VELOCITY  DATA:  PETRIE 
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TABLE  6 — Continued 
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Fig.  9.  Calculated  rotation  curves:  Petrie.  The  results 
are  listed  in  table  6.  The  value  t (mid-eclipse )  =  1.827 
days  was  used  to  calculate  6,  the  longitude  from  mid¬ 
eclipse.  The  curves  are  based  on  different  values  of  V  : 

(a)  V  =160  km/s;  (b)  Ve  =  180  km/s;  (c)  V e  =  200  km/s? 

Open  circles  represent  observations  from  the  DAO;  crosses 
indicate  a  DDO  plate. 
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in  the  data,  and  the  effect  seems  to  be  slightly  asymmetric. 
Petrie's  analysis  does  not  include  the  effects  of  limb-dar¬ 
kening,  and  assumes  that  the  axis  of  rotation  is  perpendicu¬ 
lar  to  the  orbit  plane.  The  value  of  is  about  180  km/s 
judging  solely  from  the  maximum  and  minimum  calculated  O-C 
values.  Possibly  due  to  other  previously  mentioned  effects, 
the  calculated  and  observed  curves  are  offset  from  one  an¬ 
other  . 


Hosokawa-Kopal  Analysis 

General  Theory 

Kopal  (1942;  1945)  rederived  formulae  for  the  ro¬ 
tational  effect  using  a  more  general  analysis,  and  extended 
it  to  arbitrarily  limb-darkened  spherical  stars.  His  solu¬ 
tions  were  based  on  the  assumption  that  the  star  rotates  as 
a  rigid  body  about  an  axis  coincident  with  the  pole  of  the 
orbit.  However  Kopal  notes  that  the  two  axes  need  not  nec¬ 
essarily  be  coincident,  and  Hosokawa  (1953)  used  Kopal ' s 
newly  developed  associated  a-functions  to  obtain  formulae 
for  binary  systems  of  spherical  stars  whose  axes  of  rotation 
are  inclined  from  the  perpendicular  to  the  orbital  plane. 

For  spherical  stars  rotating  as  rigid  bodies  about 
a  fixed  axis,  the  observed  radial  velocity  of  a  star  has  two 
parts:  V,  due  to  space  and  orbital  motion;  and  6V  due  to  the 
star's  axial  rotation.  The  amount  6V  is  found  by  multiplying 
the  line-of-sight  radial  rotational  velocity  v(x,y,z)  at  any 
point  on  the  stellar  surface  by  the  corresponding  element  of 
light  d£,  integrating  over  the  visible  portion  of  the  eclipsed 
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disc,  and  dividing  by  the  total  light  from  the  star.  The 
origin  of  the  rectangular  co-ordinate  system  is  the  centre 
of  mass  of  the  primary  star,  i.e.  the  eclipsed  star.  The 
X-axis  is  constantly  in  the  direction  of  the  centre  of  the 
secondary  star  projected  onto  the  xy-plane  (perpendicular  to 
the  line-of-sight) ,  and  the  Z-axis  is  along  the  line-of-sight . 
Then 


6V 


f  fV  d£ 

77dT 


(4.16) 


For  a  rigid  rotator  with  constant  angular  velocity  go,  the 
radial  rotational  velocity  V  at  any  point  is  proportional  to 
the  distance  of  the  point  from  the  axis  of  rotation  projected 
onto  the  celestial  sphere  (xy-plane) .  Figure  9  shows  the 
configuration  and  defines  the  various  parameters.  Then 


V  =  cosina{xcos  (n+3)  +  ysin(n+3))  (4.17) 


For  most  stars,  a  linear  law  of  limb-darkening  has  been  found 
to  approximate  conditions  with  a  high  degree  of  accuracy; 
so  the  distribution  of  surface  brightness  is: 

J  =  J  (1  -  u  +  ucosy)  (4.18) 

o 

where  u  is  the  co-efficient  of  limb-darkening  of  the  eclipsed 
star;  y  is  the  angle  of  foreshortening,  i.e.  the  angle  be¬ 
tween  a  surface  normal  and  the  line-of-sight  or  Z-axis;  and 

J  is  the  brightness  at  the  centre  or  the  intensity  of  rad- 
o 

iation  emerging  normal  to  the  surface. 
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Orbital  plane 


Fig.  10.  Configuration  and  notations  for  Hosokawa's  analy¬ 
sis  (from  Hosokawa  1953,  fig.  1)  with  the  following  explana¬ 
tion:  (diagram  shows  the  situation  just  after  mid-eclipse) 

(x,y):  rectangular  co-ordinates  in  plane  perpendicular 

to  the  line-of-sight ,  origin  at  the  centre  of 
the  eclipsed  disc,  x-axis  in  the  direction  of 
the  projected  centre  of  the  companion 


i:  the  inclination  of  the  orbital  plane  to  the 

celestial  sphere 


a:  the  inclination  of  the  axis  of  rotation  to  the 

line-of-sight 

n:  the  angle  between  the  y-axis  and  the  pole  of 

the  orbit,  projected  to  the  celestial  sphere 

3:  the  angle  between  the  projected  axis  of  rotation 

and  the  projected  pole  of  the  orbit 

0:  the  phase  angle  of  the  companion  in  the  relative 

orbit  (angle  ZOC) 

6:  the  projected  separation  of  the  centres  of  the 

stars 


Directions  of  rotation  and  revolution  are  the  same 
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The  limits  of  the  integration  are  briefly  determined  as  fol¬ 
lows.  The  primary's  projected  disc  is  given  by: 

2  2  2 

x  +  y  =  R±  (4.19) 

while  that  of  the  secondary  is: 

(6-x)^  +  y^  =  or  y^  =  -  (6-x)^  (4.20) 

where  =  the  absolute  radius  of  the  primary 

R^  =  the  absolute  radius  of  the  secondary 
6  =  the  projected  distance  between  the  centres  of  the 
two  stars. 

During  partial  eclipse,  the  two  discs  overlap  and  the  x- 
co-ordinate  of  the  common  chord  joining  the  two  points  of 

intersection  is  given  by: 

2  2  2 
R-,  -  R9  +  6 

s  =  - 2! -  (4.21) 

The  double  integral  is  split  up  into  two  parts,  each  bounded 
by  this  common  chord  and  the  appropriate  star's  arc:  QR^Q'S 
and  QR2Q'S. 


Fig.  11  .  Limits 
for  eclipse  integrals 
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Letting  denote  the  absolute  radius  of  the  eclipsed 
star,  we  see  that: 


or 


2  2 
l  =  x  +  Y  + 


2 

z 


and 


cos  y 


z 


R 


1 


(4.22) 


(4.23) 


Then  the  element  of  light  is  given  by: 


d£  =  —  {R  (1-u)  +  uz}dxdy  (4.24) 

JA  ^  X 

Combining  the  formulae  for  V  and  d£,  and  using  the  fact  that 
/ ydxdy  vanishes  when  integrated  over  any  area  symmetrical 
with  respect  to  the  Y-axis,  we  obtain: 


R^  (1-u) f  fx dxdy  +  u//xzdxdy 

6V  =  LOSinaCOS  (n+3)  — „ - n - rTr'J  j - ; - r  r  '  "j - 

R^  (1-u) / /dxdy  +  u//zdxdy 


(4.25) 


As  previously  stated,  the  integrations  extend  only  over  the 
visible  portion  of  the  star.  It  is  most  convenient  to  inte¬ 
grate  over  the  whole  disc,  and  subtract  the  integral  over 

the  eclipsed  portion.  The  numerator  vanishes  when  integrated 

2 

over  the  whole  disc,  and  the  denominator  reduces  to  ttR^  and 
2  2 

/3 ttR^ .  The  remaining  integrals  then  extend  over  the 
eclipsed  area,  and  are  of  the  form  which  Kopal  termed  his 
'associated  a-f unctions ' : 


TTR^+n+2 


c  p  m  ri  -  j 

ffx  z  dxdy 


(4.26) 
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Then  for  a  partial  eclipse  the  integral  becomes: 


m+n+2 

ttR^ 


s 


+  /r2-(6-x)2 


m  n  , 
x  z  dxdy 


(4.27) 


If  the  eclipse  is  total  (R^<  the  area  QR^Q'S  becomes 

equal  to  the  whole  area  of  the  primary  (eclipsed)  star,  while 
if  the  eclipse  is  annular  as  for  DE  Dra  (R^>  R2 )  the  sector 
QR^Q'S  vanishes  and  the  integral  becomes: 


( 6+R0  )  +AI- (6-x)  2 


„m+n+2  m  _ 
ttR-  a  = 

1  n 


m  n 

x  z  dxdy 


n - 7 

(6-R2)  -/r2~(6“x) 


(4.28) 


The  associated  a-functions  of  order  m  and  index  n 
all  turn  out  to  be  functions  of  two  independent  variables 
which  are  usually  taken  to  be  k,  the  ratio  of  radii  of  the 
two  components,  and  p,  the  geometrical  depth  of  the  eclipse: 


k 


rsmall 

rlarge 


and 


P 


6 


r 


rlarge 

small 


(4.29) 


where  r  , ,  ,  refer  to  the  fractional  radii.  By  defin- 

small, large 

ition  k  lies  between  zero  and  one,  while  p  is  normalized  to 
vary  from  +1  at  first  contact  to  -1  at  internal  tangency 
(totality  or  annular  phase) .  All  a-  functions  vanish  when 
p  =  1;  at  totality  where  R^<  R2  odd  orders  vanish,  while  zero 
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or  even  orders  reduce  to  constants.  However  as  in  DE  Dra 
where  R^>  ,  all  a-functions  remain  finite  at  the  moment  of 

internal  tangency. 

In  terms  of  the  associated  a-functions,  the  formu¬ 
la  for  the  rotation  effect  becomes  (see  equation  4.25): 

Cl-uJa1  +  ua1 

5V  =  -Rcosinacos  (ri+ 3)  - - - - - ~ -  (4.30) 

(l-u)  (l-a°)  +  uC/3  -  a°) 

Tables  of  values  which  enable  the  calculation  of  the  two 

lowest  functions  of  zero  order,  a°  and  a°,  have  been  tabulated 

o  1 

by  Zessewitch  (1939;  1940):  a°  represents  the  fractional  loss 

of  light  during  the  eclipse  of  a  uniformly  bright  (u  =  0) 

circular  disc  (equal  to  the  fractional  area  eclipsed) ; 

"V-  a°  represents  the  fractional  loss  of  light  of  a  circular 
^  1 

disc  completely  darkened  at  the  limb  (u  =  1) .  The  formulae 
for  the  loss  of  light  during  any  type  of  eclipse  can  be  set 
up  and  evaluated  in  a  closed  form.  However  the  explicit 
forms  are  so  complicated  that  they  are  of  practical  use  only 
when  their  values  are  presented  in  tabular  form.  For  prac¬ 
tical  reasons  the  loss  of  light  is  normalized  to  vary  between 
0  and  1  by  dividing  by  the  loss  of  light  at  the  moment  of 
internal  tangency.  In  Zessewitch’ s  tables  the  normalized 
values  for  arbitrary  limb-darkening  are  denoted  as  a  (k,p) . 

In  practice  all  a  's  can  be  obtained  from  the  appropriate 
combination  of  the  a(k,p)  for  a  uniform  (x  =  0)  disc  and  the 
aD(k,p)  for  a  completely  darkened  disc.  For  the  theory  and 
formulae  pertaining  to  eclipsing  binary  systems  see  Kopal 
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(1946) .  Due  to  the  normalization  the  tabulated  values  also 

depend  on  whether  the  eclipse  is  an  occultation  or  a  transit. 

2 

For  a  uniform  disc  the  two  are  related  by  k  ,  but  for  the 
darkened  case  the  relationship  is  more  complicated,  and 
separate  tables  of  aD(k,p)  must  be  used  depending  on  whether 
an  occultation  (function  a’)  or  a  transit  (function  a")  is 
observed . 

The  associated  a-functions  are  a  prerequisite  for 
interpreting  the  rotational  effect  in  close  binary  systems, 
for  predicting  non-radial  light  or  velocity  changes,  and  for 
interpreting  line  profiles  of  distorted  stars  during  and  out 
of  eclipse,  as  well  as  many  other  phenomena  observed  in  close 
binary  systems.  Kopal  investigated  these  functions  in  sev¬ 
eral  papers  (see  e.g.  Kopal  1942b;  1945)  and  in  1947  published 
a  detailed  study  of  the  theory  and  tabulated  values  of  this 
family  of  functions.  He  found  that  the  a-functions  of  order 

m  =  0,1,2, 3  could  be  expressed  as  simple  combinations  of 

x*  6 

certain  'I-integrals'  factored  by  powers  of  2/  and  / 

rl  r2 

or  S/  .  The  I-integrals  involve  only  one  independent  vari- 
r  2 

/  ^  \ 

able  y  =  _  ,  and  so  are  much  easier  to  tabulate. 

r2 

For  each  set  of  indices  I™  values  are  calculated  for  the 

p  /  yr 

range  of  y,  once  for  partial  phases  and  once  for  annular 

phases.  From  the  tabulated  values  of  1°  the  values  of  a1 
^  i  ,  o  o 

at  each  phase  are  found,  while  values  of  1° 

i/i 

appropriate  values  of  a*. 


yield  the 
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Calculations  of  a~Functions  for  DE  Draconis 

Because  of  the  large  number  of  eclipse  data  points 

and  the  obvious  groupings  at  various  phases,  it  was  decided 

to  reduce  the  53  observations  to  a  more  easily  managed  set 

of  15  normal  points  (simple  averages)  as  shown  in  figs.  12,13 

The  calculations  of  the  required  a°,  a° ,  a1,  a1 

o  i  o  1 

were  performed  for  each  normal  phase  point,  in  the  following 
steps,  with  values  of  k,  rgarge'  rsman'  an<^  u  kein<3  taken 
from  the  photometric  results  of  Fuertig  and  Mueninger  (1976). 

1.  For  each  point,  the  primary's  phase  in  days  was  converted 
to  true  anomaly  in  the  relative  orbit,  v,  using  the  spectro¬ 
scopic  value  of  eccentricity  and  standard  formulae  (see  e.g. 
Aitken  1964;  Schlesinger  1908) 

2.  From  this  calculated  value  of  v,  the  true  anomaly  of  the 
secondary  in  the  relative  orbit  reckoned  from  mid-primary 
minimum,  0,  was  determined  from  the  relation  (recall  that 

co  is  the  longitude  of  periastron)  : 

cos0  =  sin(v+co)  (4.31) 

3.  The  value  of  6,  the  projected  separation  of  the  centres 
of  the  two  stars,  is  given  by: 

o  o 

(4.32) 


.2  .  20  .  2.  2. 
6  =  sm  0sin  i  +  cos  i 


where  i  is  the  inclination  of  the  pole  of  the  orbit  from  the 


line-of -sight . 

4.  The  value  of  the  parameter  p  [  = 


(6-r*> 


]  was  calculated, 


and  Zessewitch's  tables  (no.  50,  1940)  yielded  the  value 
a(k,p)  for  a  uniformly  bright  disc  from  which  was  found 
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PHASE  (days) 


Fig.  12.  Eclipse  data  points.  Crosses  represent  DDO 
plates;  open  circles,  DAO  observations. 


■ 

Fig.  13.  Normal  (averaged)  points  from  eclipse  points 
of  fig.  12. 


77 


o 

c\i 


CO 


00 


U 

I 

O 


o 

OJ 


£ 

os: 


PHASE  (days) 


78 


o 

a 

o 


a  (k,p) 


5.  For  partial  phases  i.e.  if  -1  <  p  <  +1,  Zessewitch's 

tables  (no.  45,  1939)  for  a  completely  darkened  disc  gave  a 

o  2 

value  of  a" (k,p)  and  thus  a  =  k  $ (k) a" (k,p) .  For  annular 
phases  i.e.  p  <  -1,  the  parameter  q  =  is  used  for 

2 

convenience,  and  since  a(k,p)  =  1  when  p  <  -1,  aQ  =  k  for 
points  during  annular  phase.  For  a  completely  darkened  disc 
a" (k,p)  =  1  +  A(k)X(k,q)  during  annular  phase  so  that 
a°  =  k  0(k)[l  +  A(k)X(k,q)].  The  functions  $  (k)  ,  A(k),  and 
X(k,q)  were  tabulated  by  Zessewitch  as  well  (1940) 

6.  The  evaluation  of  the  higher  order  functions,  a1  and  a1, 

o  1 

involve  the  use  of  Kopal's  tabulated  I-integrals  (1947). 

g 

This  requires  the  calculation  of  y  and  D  =  —  .  For  conven- 

r2 

ience  Kopal  introduces  a  new  variable,  a,  which  is  an  inde¬ 
pendent  angular  variable  and  a  single-valued  function  of  y. 
For  partial  eclipses  (r^t^  >  6  >  r^-^),  a  =  cos  ^y .  For 

annular  phases  (6  <  Ir^-^l),  a  =  2sin  -  .  Then  the 

i-y 


functions  are  given  by: 

a1  =  2k3DI°  „ (a) 


i 


4  3/2 

and  a1  =  2k  D  1°  (a) 

l  ni 


(4.33) 


The  values  of  the  I's  are  given  at  intervals  of  5°  in  a,  so 
to  aid  in  more  precise  interpolation  than  simple  linear,  the 
modified  second  tabular  difference  M"  is  given  as  well.  The 
first  four  terms  of  the  Everett  formula  were  then  used  to 
find  the  I-values  for  intermediate  angles: 
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I  =  (l-n)I  +  nl.  +  M " E "  +  Ml'E?  (0  <  n  <  1) 
n  o  1  oo  11  —  — 

The  values  of  the  co-efficients  E"  and  E''  were  taken  from 

o  1 

Thompson  (1921) . 


The  Final  Formula 

From  the  photometry  of  Fuertig  and  Mueninger,  the 
limb-darkening  co-efficient  u  is  taken  to  be  0.5.  On  sub¬ 
stituting  this  into  the  rotation  effect  equation  4.30,  it 
becomes : 

a1  +  cl 

6V  =  -Riosinacos  (n+3)  -r - - - - -  (4.34) 

O  0  0 

■**  -  a  -  a 

5  o  i 


Also  from  the 

cost 


spherical 


sinG 

6 


triangle 

and 


relations : 

cosicosG 
sing  =  - t - 


(4.35) 


Then  setting  the  linear  equatorial  velocity  Vg  =  wR  the  above 
expression  4.34  becomes: 


6V  =  AF cos 0  -  BFsinG 


(4.36) 


where  A  =  V  sinasinftcosi  and 
e 

dependent  of  phase,  while 


B  =  V  sinacosg 
e 


are  both  in- 


F 


1 

6 


a 1  +  a  1 
o  1 


(4.37) 


is  dependent  on  phase.  The  values  of  the  factors  involved 
in  F,  and  the  results  for  each  phase  angle  are  presented  in 
Table  7. 

In  order  to  fit  the  observed  velocities  with  some 


ROTATION  EFFECT:  HOSOKAWA 
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calculated  points,  various  values  of  3  were  substituted  and 
an  average  value  of  Vasina  was  found  using  the  observed 
<5Vobs*  Because  the  0-C  curve  is  only  slightly  skew,  and  the 
velocity  curve  itself  is  quite  symmetric,  the  axis  of  rota¬ 
tion  cannot  be  radically  inclined  from  the  pole  of  the  orbit. 
Thus  we  are  restricted  to  small  values  of  $,  and  values  of 
a  close  to  i  [-  87°].  The  average  value  Vesina  =  188  km/s 

was  used  to  calculate  the  expected  6V  ,  for  each  normal 

caic 

phase  point  for  3=0°,  +5°,  ±10°.  From  the  calculated 

curves  in  Figure  14  compared  with  the  observed  points  the 

best-fit  value  of  3  to  the  maximum  and  minimum  O-C  is  a  small 

positive  quantity  of  about  +5°.  As  with  Petrie’s  analysis, 

the  calculated  curves  are  off-set  from  the  observed  normal 

points,  partly  due  to  some  other  effects  not  included  in 

this  analysis,  e.g.  gravity-darkening. 

The  value  of  a  should  not  reasonably  be  less  than 

about  70°,  so  from  V  sina  =  188  km/s,  the  value  of  V  could 

e  e 

be  between  about  177  and  188  km/s.  By  Petrie’s  previous 
method,  his  value  of  about  180  km/s  falls  right  in  the  middle 
of  the  range  expected  from  Hosokawa’s  analysis.  If  a  is 
taken  to  be  the  inclination  of  the  orbit,  -  87°,  then 
Ve  =  188  km/s  which  is  slightly  higher  than  Petrie's  value. 

Comments 

The  values  from  both  analyses  above  fall  within 
the  range  expected  for  a  single  B9V  star.  However  if  syn¬ 
chronism  is  assumed,  then  the  should  only  be  about  30  km/s 
assuming  the  radius  to  be  3.1  solar  radii.  Obviously  the 


# 

, 
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Fig.  14.  Calculated  rotation  curves:  Hosokawa.  The  data 
are  listed  in  table  7.  The  four  curves  are  based  on  the 
same  value  of  Vesina  =  188  km/s,  but  different  values  of  3: 

(a)  3  =  -10°;  (b)  3=0°;  (c)  3  =  +5°;  (d)  3  =  +10°. 

Open  circles  represent  the  normal  (averaged)  points  from 
the  observational  data. 
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system  does  not  rotate  in  synchronism,  and  so  the  stars 

should  rotate  like  single  stars,  thus  supporting  the  value 

for  V  obtained  for  DE  Dra. 
e 


CHAPTER  V 


CONCLUDING  REMARKS 

Thus  far,  this  paper  has  been  concerned  with  the 
physical  parameters  of  the  spectroscopic  binary  system  DE 
Dra  which  can  be  obtained  from  the  radial  velocity  measure¬ 
ments  of  spectrograms,  as  well  as  information  on  the  indi¬ 
vidual  components  deduced  by  combining  the  spectroscopic 
results  with  some  preliminary  photometric  data.  Many  of  the 
results  had  to  be  calculated  using  assumed  values  for  a  nor¬ 
mal  B9V  star,  and  so  are  only  approximate.  However  they  do 
serve  to  support  some  general  theories  concerning  this  class 
of  early-type  star  and  in  that  respect  even  the  rough  calcu¬ 
lations  serve  a  useful  purpose. 

As  with  all  scientific  research,  though,  there  are 
many  areas  which  require  further  investigation  in  order  to 
build  up  a  complete  and  firmer  knowledge  of  this  system.  In 
the  area  of  spectroscopic  studies  the  fairly  detailed  analy¬ 
sis  presented  here  has  dealt  only  with  the  velocity  curve  of 
the  primary,  although  the  number  and  quality  of  the  plates 
has  established  this  very  reliably.  High  dispersion  spectro¬ 
grams  in  the  red  or  infrared  taken  near  the  extremes  of  the 
velocity  curve  (that  is,  at  quadratures)  may  allow  the  sec¬ 
ondary's  spectrum  to  be  seen  and  thus  allow  its  velocity 
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curve  to  be  calculated.  With  velocity  curves  from  both  com¬ 
ponents,  certain  orbital  elements  (a,  i)  can  be  calculated 
and  compared  with  those  from  photometric  data. 

Some  recent  higher  dispersion  plates  taken  during 
eclipse  before  and  after  conjunction  do  show  evidence  of 
double  lines.  However  this  effect  is  thought  to  be  due  to 
the  differential  rotational  velocity  of  the  visible  limbs: 
the  lines  are  double  because  although  part  of  the  centre  of 
the  star  is  obscured,  the  limbs  which  are  approaching  or  re¬ 
ceding  are  still  visible.  A  measurement  of  the  splitting 
could  serve  as  a  means  of  calculating  the  absolute  rotational 
velocity  since  there  would  be  no  overall  blending  from  the 
whole  star.  The  line  shifts  are  due  to  each  limb  alone. 

However,  higher  resolution  is  needed  before  any  reliable 
measurements  of  the  separation  can  be  made.  It  appears  as 
well  that  this  may  be  the  first  system  in  which  this  aspect 
of  the  separation  of  different  rotation  velocities  at  the 
limbs  has  been  seen  and  recognized. 

The  photometric  results  for  DE  Dra  are  preliminary 
only.  A  more  detailed  analysis  with  data  covering  more  than 
one  eclipse  would  yield  better  values  of  the  inclination, 
the  radii,  and  the  limb-darkening.  In  particular,  if  the 
light  curve  is  observed  in  more  than  one  effective  wavelength 
the  co-efficient  of  limb-darkening  can  be  more  accurately  de¬ 
termined.  This  results  from  the  fact  that  darkening  not  only 
decreases  the  brightness  from  centre  to  limb,  but  also 
causes  the  limb  to  be  redder  than  the  centre.  Thus  during 
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an  eclipse  the  colour  index  of  the  darkened  star  should  show 
a  gradual  change.  It  should  be  noted  here  that  any  reflec¬ 
tion  effect  in  the  system  (which  has  not  been  considered  in 
this  paper)  will  tend  to  exaggerate  the  limb-darkening. 

Due  to  the  great  difference  in  magnitudes  and  sizes 
of  the  two  components,  it  is  not  likely  that  reflection  or 
tidal  distortion  are  important  effects  for  the  primary.  How¬ 
ever  the  following  facts  could  be  kept  in  mind.  The  reflec¬ 
tion  effect  shifts  the  apparent  centre  of  light  of  each  com¬ 
ponent  in  the  direction  of  the  centre  of  mass  of  the  system. 
The  effect  vanishes  at  conjunction  and  is  a  maximum  at  quad¬ 
ratures.  If  not  corrected  for,  this  could  diminish  the  true 
amplitudes,  K,  of  the  velocity  curves  and  lead  to  an  under¬ 
estimate  of  the  true  masses  and  dimensions  of  the  stars. 

Tidal  distortion  also  affects  K  and  reinforces  the  reflection 
effect  if  i  <  63?4  but  counteracts  it  if  i  >  63?4.  Also 
ellipticity  (of  stellar  shape)  and  reflection  will  cause  the 
radial  velocity  curves  to  be  asymmetric,  producing  a  spurious 
orbital  eccentricity. 

An  error  in  the  limb-darkening  will  affect  the 
analysis  of  the  rotation  effect.  As  well,  the  discrepancy 
between  the  calculated  curves  and  the  normal  points  show 
that  for  a  complete  picture,  effects  other  than  just  limb- 
darkening  must  be  included.  These  encompass  gravity-dark¬ 
ening,  since  to  every  degree  of  limb-darkening  there  corres¬ 
ponds  a  certain  amount  of  gravity-darkening.  Centrifugal  and 
tidal  forces  cause  stars  to  deviate  from  spheres  (an  assump- 
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tion  in  Chapter  IV)  and  have  non-uniform  surface  brightness. 
This  is  due  to  the  fact  that  limb-darkening  makes  brightest 
the  parts  of  the  surface  nearest  the  observer,  while  gravity- 
darkening  makes  brightest  those  nearest  the  star's  centre. 

In  general,  on  a  distorted  star  the  point  closest  to  the  ob¬ 
server  is  not  the  geometric  centre  of  its  apparent  disc. 

In  connection  with  this,  gravity-darkening  at  the  tidal 
bulge,  which  usually  lags  behind,  can  cause  the  times  of 
spectroscopic  conjunction  to  precede  mid-eclipse  times,  of 
which  there  is  some  evidence  in  the  observations  of  DE  Dra. 

The  variation  of  angular  velocity  from  pole  to  equator  (in 
Chapter  IV  co  was  assumed  constant)  should  hlso  be  included 
in  the  rotation  effect  analysis.  The  theory  and  equations 
needed  for  all  these  are  available,  and  it  remains  to  combine 
and  apply  them  to  DE  Dra  to  obtain  a  more  precise  picture  of 
its  rotation. 

In  addition  to  velocity  measurements,  analyses  of 
observations  of  stellar  radiation  with  high-dispersion  spec¬ 
trograms  can  lead  to  a  knowledge  of  the  conditions  (pressure, 
temperature,  chemical  composition)  in  the  stellar  atmosphere. 
The  available  spectrograms  also  have  an  intensity  calibra¬ 
tion  spectrum,  so  measurements  can  be  made  which  would  give 
line  profiles  and  equivalent  widths.  The  spectrum  is 
scanned  with  a  microphotometer  to  obtain  a  density  versus 
wavelength  tracing.  This  is  converted  into  an  intensity 
scale  through  the  intensity  calibration  spectrum  on  the  plate. 
The  observations  are  expressed  relative  to  an  adopted  datum 
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line  called  the  continuum  of  the  star,  which  for  early-type 
stars  with  few  strong  lines  and  many  weak,  shallow  ones,  is 
quite  easy  to  define.  This  represents  the  intensity  distri¬ 
bution  of  a  stellar  atmosphere  at  the  same  temperature  and 
pressure  as  the  real  one,  but  without  any  line  absorption 
just  absorption  in  the  continuum.  The  continua  observations 
can  be  compared  directly  to  continuous  fluxes  predicted  by 

models  to  determine  a  model  T  It  should  be  noted  that 

ef  f 

any  analysis  must  use  rectified  intensity  tracings,  i.e. 
corrected  for  distortions  due  to  interstellar  matter,  Earth's 
atmosphere,  and  instrumental  broadening.  Line  profiles  are 
drawn  from  the  rectified  tracings  and  the  area  of  the  profile 
determined,  which  will  yield  the  equivalent  width  of  the 
line.  The  equivalent  width  of  a  line  increases  with  in¬ 
creasing  element  abundance  and  its  distribution  among  the 
various  states  of  ionisation  and  excitation.  The  shape  of 
the  curve  of  growth  changes  with  atmospheric  temperature 
and  pressure.  From  the  spectrum  then  a  set  of  parameters 
for  the  observed  curve  of  growth  can  be  converted  through 
well-investigated  formulae  into  information  about  atmospheric 
temperature,  pressure,  state  of  atmospheric  motion,  and 
element  abundance. 

Alternately,  the  observed  line  profiles  and  equi¬ 
valent  widths  can  be  compared  with  those  predicted  by  model 
atmospheres  which  are  characterised  by  the  parameters  T 
surface  gravity  g,  and  chemical  composition.  In  early-type 
stars  the  strengths  of  the  Balmer  lines  are  both  temperature 
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and  pressure  dependent:  the  outer  wing  shape  is  determined 
by  the  density  of  the  ions,  the  inner  portion  reflects  the 
temperature  structure.  Thus  equivalent  widths  (element  abun¬ 
dances)  can  be  estimated  by  comparing  the  observed  curve 
with  a  set  of  theoretical  profiles  having  the  same  shape  but 
different  areas.  A  value  for  Teff  is  obtained  by  comparing 
the  observed  Balmer  jump  with  that  predicted  by  models.  Also 
if  is  correct,  the  relative  strengths  of  the  lines  for 

different  stages  of  ionisation  of  one  element  should  yield 
the  same  element  abundance.  A  knowledge  of  the  element  abun¬ 
dances  in  the  atmospheres  of  normal,  main  sequence,  early- 

type  stars  is  important  because  they  are  among  the  youngest 

6  8 

objects  in  the  HR  diagram  (about  10  to  10  years  old) . 

Their  surface  compositions  can  thus  provide  information  about 

the  present  chemical  composition  of  the  interstellar  matter 

they  were  formed  from,  and  a  comparison  with  older  objects 

9 

(the  Sun,  for  example,  is  about  4.5x10  years  old)  could 
allow  a  calculation  of  the  rate  of  enhancement  of  interstellar 
matter  with  heavier  material. 

As  a  side  effect  of  obtaining  line  profiles  and 
widths  from  the  spectra  of  DE  Dra,  an  independent  determin¬ 
ation  can  be  made  of  the  rotational  velocity,  which  would 
prove  interesting.  Axial  rotation  shifts  the  cores  of  the 
spectral  lines,  but  also  widens  the  line  profiles.  The  pro¬ 
files  are  asymmetric  while  the  star  is  being  eclipsed,  and 
an  analysis  of  the  extent  and  form  of  the  line  asymmetry  can 
give  a  value  for  the  rotational  velocity  (and  also  the  geo- 


' 


91 


metry  of  the  eclipse  itself) . 

Although  the  spectroscopic  analysis  of  the  primary 
component  of  DE  Dra  is  fairly  complete,  it  is  evident  that 
much  more  information  can  be  obtained,  not  only  from  further 
investigations  on  present  spectral  data  but  also  from  inves¬ 
tigations  in  other  observational  and  theoretical  areas.  It 
is  hoped  that  such  information  as  presented  here  will  aid 
and  stimulate  further  research  on  this  system. 
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APPENDIX  A 


OBSERVATIONAL  DATA 

The  observations  presented  here  include  both  non¬ 
eclipse  and  eclipse  data  from  the  DDO  and  the  DAO.  Obser¬ 
vations  numbered  1  to  19  inclusive  are  from  DDO  plates,  of 
which  numbers  18  and  19  are  eclipse  data.  The  remaining 
observations  from  numbers  20  to  105  are  from  DAO  plates, 
and  of  these  the  last  51,  from  numbers  55  to  105  inclusive, 
are  eclipse  data. 

The  column  headings  denote  the  following 


quantities : 

V: 

the  observed  heliocentric  radial  velocity 
in  km/s 

TIME: 

the  phase  of  the  observation  in  days 

O-C : 

the  observed  velocity  minus  the  calcu¬ 
lated  velocity  in  km/s 

MEAN  A: 

2  Ti¬ 
the  mean  anomaly  M,  in  radians,  =  — (t-T) 

ECC  A: 

the  eccentric  anomaly  E,  inradians,  de¬ 
termined  from  the  formula  M  =  E-esinE 

TRUE  A: 

the  true  anomaly  v,  in  radians,  calculated 

from  ,  ,  t-, 

.  v  1+e  ,  E 

tan  -Tc  =  ■= -  tan 

2  1-e  2 

U: 

the  argument  of  the  latitude  u  =  v+cof  in  rad 

C: 

the  velocity  in  km/s  calculated  from  the 
orbital  elements,  =  y  +  Kecosoo  +  Kcosu 
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APPENDIX  B 


CALCULATION  OF  PHYSICAL  PARAMETERS 
Colour  and  magnitude 

The  description  of  the  uvbyB  system  given  by  Crawford 

(1975)  contains  the  following  definitions  of  the  parameters: 

E(C;l)  =  0 . 20E  (b-y ) 

E (u-b)  =  1.5E(b-y) 

m^  =  (v-b)  -  (b-y) 

c1  =  (u-v)  -  (v-b) 

(u-b)  =  cg  +  2m^  +  2  (b-y) 

(u-b)  =  1.13  (u-b)  -  1.70  (b-y)  -  0.186 

o 

The  photometry  of  Crawford  et  al  (1973)  resulted  in  values 
of  the  indices  for  DE  Dra  of: 

8  =  2.824 
b-y  =  -0.017 
m^  =  0.123 
c  =  0.873 

These  values  are  substituted  into  the  above  equations  to  yield 

(u-b)  =  1.085 

(u-b)  =  1.069 

o 

Then  E  (u-b)  =  0.016  [=  (u-b)  -  (u-b)  ] 

o 

and  E(c^)  =  0.002 
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The  subscript  ’o’  and  the  letter  E(  )  refer  to  the  unreddened 

index  and  the  colour  excess  respectively.  Thus  the  unreddened 

parameter  cq  can  be  found.  This  allows  us  to  calculate  the 

absolute  visual  magnitude  using  the  recent  calibration  of 

the  uvbyB  system  performed  by  Crawford  (1978).  He  gives  a 

table  relating  the  c  value  to  a  value  of  3 (ZAMS)  (i.e.  the 

o 

strength  of  the  H3  line  for  unevolved  B-type  stars) ,  which 

is  then  related  to  a  value  for  the  M  (ZAMS) .  However,  due 

v 

to  evolutionary  effects,  the  value  of  3  (observed)  will  not 

be  the  same  as  the  tabulated  3 (ZAMS)  for  that  value  of  c  . 

o 

Therefore  a  correction  must  be  made  to  the  M  (ZAMS)  to  ob- 

v 

tain  the  correct  M  .  The  above  calculations  for  DE  Dra  pro- 

v 

cede  as  follows: 

cq  =  c^  -  E(c^)  =  0.871  (from  previous  page) 

3 (ZAMS)  =  2.857  (from  Crawford  1978) 

M  (ZAMS)  =  1.08 
v 

A3  =  3  (ZAMS)  -  3  (observed)  =  0.033 

AM  =  10A3  =  0.33 
v 

Then 

M  =  M  (ZAMS)  -  AM  =0.75 
v  v  v 

This  value  agrees  well  with  that  expected  of  a  B9V  star. 

Crawford  also  gave  a  relation  between  the  uvby3  system  and 

the  UBV  system: 

(U-B)  =  -0.94  +  0.86c 

o  o 


■ 
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This  gives  a  value  for  DE  Dra  of  (U-B)  =  -0.19  which  also 

o 

agrees  well  with  the  value  expected  for  a  B9V  star  in  the 
MK  system. 


Mass 

Different  values  for  the  mass  have  been  calculated 
from  the  various  formulae  given  by  five  authors.  The  final 
value  adopted  in  chapter  III  is  an  average  of  these. 

1.  Allen  (1973) 


(for  the  value  of  M,  see  the  sections  Bolometric 

bol 

correction  and  Bolometric  magnitude) 

M  _  =  4.75  -  2 . 51og  k 
bol  L 


© 


log  =  3 . 451og  ~ 
©  © 


Then 


M  =  3.38M 


© 


2.  Batten  (1968) 


M^  =  -5.5591og(M)  +  2.485 


M  =  2.05M 


© 


3.  Petrie  (1950) 


log  (M )  =  0.548  -  0.144M  +  0.002M' 


M  =  M,  ,  +  2  log  e  =  0.8  64 
c  bo1  5T07 


M  =  2.66M 

© 

4.  Harris,  Strand,  and  Worley  (1963) 

M,  ,  =  4.6  -  10.01og(  ) 
bol 

M  =  2.76M 

© 


for  M 


bol 


+  7.5 
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5.  Kopal  (1959) 

log  (M)  =  0.45  -  0.143M,  ,  for  M  >  2m 

bol  © 

M  =  2.65m 

© 

The  average  value  of  the  five  values  above  is  then 

M  =  2.70M 
av  © 


Distance 

The  value  of  V,  the  apparent  visual  magnitude  uncor¬ 
rected  for  reddening,  was  published  by  Crawford  et  al  (1973) 
for  DE  Dra  as  V  =  +5.72.  The  standard  relationship  for  inter¬ 
stellar  absorption  was  used,  as  well  as  the  values  from  the 
photometry  for  the  colour  excesses. 

E(b-y)  =  0 . 7  4E  (B-V ) 

E  (u-b)  =  1.5E(b-y) 


Then 


(Crawford  1975) 


A  =  3 . 3E  (B-V )  =  T^T  E  (u-b) 

v  .741.5 

=  +0.048 


V  =  V  -  0 . 048 
o 

=  +5.67 


V  -  M  =  51og(d)  -  5  =  +4.92 
o  v  ^ 


d  =  96.4  parsecs 


Effective  temperature 

The  effective  temperature  of  11,300°K  is  an  aver¬ 
age  value  of  three  quantities  obtained  from  formulae  given 

5040 

by  Heintze  (1973)  relating  (U-B)  to  0  =  = - . 

°  e“  eff 


. 
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Bolometric  correction 

The  bolometric  correction  was  calculated  using  the 
formula  in  Allen  (1973) : 


-B. 


c . 


=  -42.54  + 


logTeff  + 


29000 

T 

eff 


B.C.  =  -0?56 


Bolometric  magnitude 


From  the  above  B.C.  the  value  of  M,  n  is: 

bol 

M  =  M  +  B.C. 
bol  v 


=  +0* I 2?19 


Radius 

An  average  value  of  the  radius  was  calculated  from 
two  formulae  relating  ^  and  R: 

1.  Allen  (1973) 

log  |  =  |  [  42.36  -  lOlogT  -  M^] 

© 


R  =  2.13  R 

© 

2.  Kopal  (1959) 

0.4  5  -  0 . 14  3M,  _ 

bol 

R  =  2. 32  R 

© 


1.571ogR  -  0.15 


for  M>2M 

© 


This  gives  a  value  for  R  of  the  primary  of: 

R  =  2.2  R  =  R-i 
av  ©  1 

and  from  the  photometric  value  k  =  0.35: 

R0  =  0.8  R 
2  © 

However  because  of  the  great  uncertainty  in  the  value  of 

M,  ,  it  was  decided  to  use  the  spectroscopic  value  of  the 
bol 
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mass  function  f  (M)  ,  the  calculated  mass  M  ,  and  the  photo¬ 
metric  value  of  i  to  calculate  values  for  a^, M  ,  and 

R0 .  From  =  2.7  M  and  i  -  87°,  r,  =  0.157,  r~  =  0.055: 

2  1  ©  '1  A 

3  .  3. 

MLsm  i 

f  (M)  =  — - ^  =  0.074 

<m1+m2) 


Mo  =  1.0  M 


© 


a^sini 


=  3.73x10 


= 


3.74x10  km 

M1  6 
— —  a,  =  9.99x10  km 
m2  l 


a  =  a^  +  a2  =  13.7x10  km 


R.,  =  ar,  =  3.1  R 
11  © 


R„  =  ar~  =  1.1  R 
2  2  © 


